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ABSTRACT 

The  analysis  oi  the  data  for  the  International  Seismic  Month 
is  nearing  completion.  The  effects  of  crust  and  upper- 
mantle  structure  on  the  short-period  signals  received  at 
LASA  and  NORSAR  are  des  •“•'bed.  Further  work  on  inter¬ 
fering  Rayleigh  waves  and  or  the  phase  pP  is  included,  to¬ 
gether  with  several  studies  «n  global  seismology.  Our  de¬ 
velopment  of  a  seismic  computer  terminal  system  for  the 
ARPA  Network  is  progressing  satisfactorily. 
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SUMMARY 


This  is  the  twentieth  Semiannual  Technical  Summary  report  of  M.I.T.  Lincoln  Laboratory's 
Seismic  Discrimination  Group.  The  research  activities  of  this  Group  are  directed  toward  the  so¬ 
lution  of  the  seismological  problems  associated  with  the  detection,  location  and  identification  of 
earthquakes  and  nuclear  explosions.  We  are  particularly  concerned  with  the  analysis  of  global 
seismic  information,  utilizing  modern  interactive  computer-processing  techniques,  in  order  to  de¬ 
velop  the  technological  background  necessary  for  the  negotiation  of  a  comprehensive  test-ban  treaty. 

A  major  effort  of  the  Group  has  been  to  study  the  detection  and  location  capability  of  the 
existing  network  of  seismic  stations  and  arrays.  The  period  20  February  to  19  March  1972  was 
designated  cs  the  International  Seisniic  Month  and,  with  the  excellent  cooperation  of  seismologists 
in  many  countries,  the  available  seismic  data  for  this  period  have  been  analyzed  in  detail.  The 
first  phase  of  this  project  consisted  of  the  preparation  of  an  epicenter  list  based  on  the  short- 
period  arrivals  at  selected  WWSSN  stations,  Canadian  network  stations,  and  a  number  of  arrays 
including  IJ^SA  and  NORSAR.  This  phase  is  now  completed,  and  the  final  epicenter  list  contains 
996  entries.  Special  attention  has  been  given  to  the  precision  of  the  estimates  of  epicentral  loca¬ 
tion,  focal  depth  and  magnitude,  and  the  events  have  been  categorized  by  the  reliability  of  these 
estimates.  The  second  phase  of  the  study  concerns  the  association  of  long-period  arrivals  with 
the  listed  events,  and  the  assignment  of  surface  wave  magnitudes.  The  problem  of  associating 
the  long-  and  short-period  data  has  raised  a  number  of  interesting  questions  that  are  important 
for  seismic  discrimination.  This  phase  will  be  completed  shoT-tly.  The  final  report  on  this 
study,  including  basic  statistical  analyses  of  the  data,  will  be  available  in  the  spring  of  1974. 

We  continue  to  have  a  strong  interest  in  research  concerning  the  performance  of  seismic 
arrays,  and  a  number  of  studies  of  this  kind  are  described.  W'e  have  developed  the  comparison 
of  NORSAR  with  LASA,  and  it  is  becoming  apparent  that  much  of  the  observed  complexity  of 
short-period  arrivals  at  NORSAR  may  be  attributed  to  crust  and  upper-mantle  structure.  Re¬ 
sults  are  presented  detailing  the  nature  of  this  complexity  for  each  subarray,  in  terms  of  spectral 
characteristics,  and  explaining  how  these  influence  the  performance  of  the  array  as  a  whole,  in 
terms  of  array  diagrams.  Some  of  these  features  are  similar  to  those  found  at  LASA. 

Research  in  seismic  discrimination  has  concentrated  on  two  subjects  -  the  interfering  Ray¬ 
leigh  wave  problem,  and  the  further  study  of  the  depth  phase  pP.  The  presence  of  multipathing 
can  complicate  the  process  of  the  separation  of  Rayleigh  waves  from  two  events,  but  results  are 
given  that  indicate  some  success  using  high- resolution  frequency-wave  number  techniques.  Com¬ 
plications  in  the  structure  of  pP  due  to  reflections  at  internal  boundaries  in  the  earth  are  de¬ 
scribed,  and  the  identification  of  pP  in  explosion  records  is  discussed. 

In  preparation  for  making  use  of  the  global  seismic  information  to  be  available  via  the  ARPA 
Network,  we  have  been  extending  our  research  into  global  seismology,  and  are  strengthening  our 
facilities  for  the  rapid  access  and  processing  of  digital  data.  Improvements  in  world-wide  seis¬ 
mic  station  corrections  are  proposed,  and  we  continue  our  study  of  inhomogeneities  in  the  earth. 
S-wave  travel-time  anomalies  show  evidence  of  these  inhomogeneities  even  more  clearly  than 
P-waves.  In  addition,  we  describe  a  study  that  is  the  first  step  in  a  series  of  investigations 
attempting  to  relate  earthquake  source  mechanisms  to  large-scale  tectonics. 

Our  development  of  a  sophisticated  interactive  computer  facil.*ty  for  utilizing  the  capabilities 
of  the  ARPA  Network  is  continuing.  We  shall  have  a  permanent  connection  into  the  network  In  the 
very  near  future,  and  we  expect  to  be  ready  to  receive  seismic  information  as  it  becomes  avail¬ 
able  oi.  the  Network. 


M.A.  Chinnery 
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S  i:  I  S  M  i  (•  1)1  sc  HI  M  I  X  AT  1  OX 


I.  [X1’1:HX  \'r!OX  VL  SHISAIIC  MOXTH 

A.  ISM  I'A'IlX’r  MSI'  C|l\HAi  r  !■:  HI  S'll  CS 

1  lit-  International  Seisinir  Alontli  flSM)  i  ;  a  nmiti-iiationa I  i  ooporat i\c‘  f\porinu;nt  wlueh  has 
been  'lisi  iisiud  in  previous  SATS  and  in  a  rei-enl  report^  Our  final  list  of  verifieil  events  for 
that  experiment  contains  entries,  of  v.  hi(  h  have  not  been  assi^neil  body  wave  ma^'nitudes. 
ri.e  most  eonitnon  reasons  for  tiiis  are- 

(1)  The'  lac'k  of  ro[K)rte(l  urnplftudes  b\'  t  ertain  stations. 

(2)  The  la»  k  of  roeorrled  P  arrivals  at  distani'es  f'roater  than  10'. 

Obviously,  the  first  of  these  also  affects  the  quality  of  the  nij^  value  even  when  we  do  Ikivc  enough 
amplitudes  to  assign  nij^  to  an  e\cut.  With  this  in  mind,  we  have  taken  steps  to  partially  correct 
the  situation.  Our  film  library  was  sear'  hefi  to  find  film  chips  for  stations  which  were  used  in 
the  hSAI  but  which  previon.siy  had  not  been  completely  r  ad  specifically  for  the  hSAl.  There  were 
20  such  stations.  A  search,  of  hSAl  data  indic-ited  that  more  than  800  arrivals  at  these  stations 
had  been  assiicrated  with  e\ents,  but  that  no  amplitude  or  period  readings  were  in  the  data  base. 
The  film  data  have  been  reviewed  and  the  missing  amplitude  and  period  data  obtained.  These 
data  will  shortly  be  integraterl  into  the  rest  of  the  hSM  information,  and  nij^  values  revi.scd  ac¬ 
cordingly.  This  is  expcctCfl  to  change  some  of  the  statistics  which  are  discussed  helow'. 

An  attempt  ha.s  been  made  to  grade  events  according  to  the  reliability  that  we  can  assign  to 
the  event  [t.ira meters.  Ihere  are  eight  grades,  namely.  A,  li,  f,  and  D  for  events  which  could 
be  locatc'd  without  restraining  the  de[)th  in  any  wav,  and  Al,  HI,  ('!,  and  1)1  for  ca.ses  when  the  data 
left  the  depth  indeterminate  iivl  it  was  nei'cssarv  tf)  arbitrarily  restrain  depth.  To  avoid  «’onfu- 
sion,  we  should  note  that  depth  phases  are  uscfl  routinely  by  our  location  pr(>gram,  and  that  a 
hyi»ocentor  obtaincrl  using  depth  phases  is  consiflered  unrestrained.  An  indeterminate  depth  indi¬ 
cates  that  the  location  pn>gram  converged  to  a  negative  depth  or  a  depth  greater  .nan  700  km.  or 
that  de[)i  w.is  so  poorly  i-onstrainefl  that  convergence  couhl  not  be  obtained. 

The  t  vent  grades  are  based  upon  six  different  parameters.  One  of  these,  N.-WKIVIH.  is  the 
number  of  indivirlual  station  ivij^  values  averagcfl  to  get  the  ISM  m^  value,  .AVKMH.  The  other 
five,  which  relate  to  location,  are 

DrCCfF  Degree  '  of  freedom  in  the  location.  This  is  the  number  of 

observations  used  minus  the  number  of  parameters  to  be 

deterrniried. 

HATH)  (loodness  of  fit  measure.  Hesiduals  are  w'cighted  inversely 

bv  their  assigned  standard  deviations,  summed,  and  divided 
by  DKfiP’  to  obtain  HATK). 

MAXA.X2  l.arg'^st  principal  semiaxis  of  our  estimated  epicenter  error 

ellipse,  (’alculatcfl  only  for  unrestrained  locations.  Units 
f  kilometers. 

DKPTIK^  i  1  length  of  our  estimated  depth  error  bar,  in  kilometers. 

(  alculated  only  for  unrestrained  locations. 

DUPTH  ISM  assigned  deptli  for  the  event,  in  kilometers. 

The  exact  quantitative  meaning  of  some  of  these  items  is  not  finalized,  but  we  have  tried  to  use 
them  nevertheless.  Some  discussion  of  them,  and  our  location  program,  will  be  found  in  our 
last  SATS  (  10  .hrne  1071).^ 
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TABLE  I-i 

LINCOLN  iSM  GRADING 

RULES 

Run  Free 

Restrained 

Grade  r  A  (287) 

MAXAX2<  100 

RATIO  <  2 

DEGF  >  4 

DEPTHQ/DEPTIU  0.5  or  DEPTHQ  <  25 
NAVEMB  >  3 

Grade  =  AI  (61) 

RATION  1 

DEGF  >  5 

NAVEMB  >  3 

Grade  -  B  (56) 

MAXAX2<  100 

RATIO  <  2 

DEGF  >  4 

DEPTHQ /DEPTHS  0.5  or  DEPTHQ^  25 
NAVEMB  =  1  or  2 

Grade  =  BI  (19) 

ration:  1 

DEGF  ^  5 

NAVEMB  =  1  or  2 

Grade  =  C  (SO) 

.MAXAX?.<  100 

RATIO  <  2 

DEGF  ^  4 

DEPTHQ/DEPTH  >  0.5 

DEPTHQ  >  25 

NAVEMB  ^  1 

Grade  =  Cl  (32) 

RATIO  >  1 

RATIO  <  2 

DEGF  >  5 

NAVEMB  1 

Grade  -  1.  .'197) 

All  tun  free  whir-h  are  not 

A.  E  or  C 

Grade  =  DI  (264) 

All  restrained  which 
are  not  A,  B,  or  C 

Table  I-i  gives  the  viefinitior.  of  grades  in  terms  of  the  above  parameters.  The  number  of 
events  in  ea^’h  grade  is  giv  'n  in  parentheses.  For  a  given  grade,  all  conditions  must  be  met. 

The  transition  from  Grafl^  A  to  D  is  fairly  clear  in  the  case  of  unrestrained  events?  that  is: 

A  -  Excellent  hypocenter  '’nd  m^^ 

B  =  Excellent  hjT>0'’enter,  bv.  perhaps  less  reliable  m^^ 

C’  r  Excellent  epic  enter,  but  peri:aps  less  reliable  depth  as  well  as  m^^ 

1)  =  All  others. 

The  transition  is  more  arbitrary  nnd  less  clear  for  AI  to  DI.  First  note  that  DEGF  ^  5  for  AI, 

BI,  and  Cl,  whereas  it  is  DEGF  5^  ^  for  A,  B,  and  G.  In  effect,  this  means  that  we  must  have 
at  least  eight  observations  used  for  location  in  all  cases. 

F’or  AI  and  BI,  we  require  RATION  1  rather  than  RATIO 2  in  the  hope  that  tighter  control 
c  1  that  parameter  will  compensate  for  not  having  DEPTHQ  or  MAXAX2.  In  case  Cl.  we  include 
events  with  i  <  RATIO <  2  in  place  of  relaxing  the  depth  control  parameters  which  aie  not  avail¬ 
able.  The  transition  in  terms  of  NAVEMB  is  uncharged  for  restrained  events.  We  should  note 
that  for  many  purposes  events  in  grades  AI.  BI,  and  Cl,  and  even  many  in  DI,  probably  have 
good  depths.  That  is,  shallow  depths  are  probably  good  within  a  few  tens  of  kilometers  and  deep 
events  do  get  listed  as  deep.  Unfortunately,  this  is  a  subjective  and  unsubstantiated  assertion, 

A  final  point  to  re-emphasize  is  that  events  were  restrained  in  depth  only  when  that  was  the  only 
way  convergence  of  the  location  program  could  be  obtained. 

Events  in  grades  D  and  DI  tend  to  be  the  lower-magnitude  events.  This  is  clear  from  Fig.  I-l. 
However  not  all  (D,  DI)  events  have  assigned  magnitudes  so  small  that  they  are  of  no  interest  for 
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discrimination.  This  creates  a  problem  unless  these  events  can  be  somehow  upgraded  or  it  can 
be  established  that  the  assigned  mag»  itude  is  very  much  in  error.  (Xir  reading  of  additional  am¬ 
plitudes  should  help  this  situation  somewhat. 

The  overall  geographical  distribution  of  all  ISM  events  is  shown  in  P'ig.  I-?.(a).  Figure  I-2(b) 
shows  the  distribution  of  events  when  (D,  DI)  events  are  excluded  from  the  population.  Table  1-2 
is  a  somewhat  more  detailed  breakdown  of  geographical  distribution  by  grade  and  seismic  region 
(Flinn  and  Engdahl^).  The  fraction  of  events  in  the  (D,  DI)  categories  for  a  reg  on  gives  a  crude 
idea  of  how  well  that  region  was  covered  by  ISM  data  and  procedures.  No  more  than  25  percent 
of  the  events  from  regions  18,  21.  24,  26.  27.  28,  ana  31  are  D  or  DI.  However  More  than 
75  percent  of  events  in  regions  3,  11,  16.  35,  36,  37,  38,  i9  and  43  are  category  L' or  DI,  lead¬ 
ing  us  to  believe  that  these  regions  require  improved  station  coverage  or  data-proc:‘Ssing  methods. 

It  is  our  belief  that  all  events  on  the  ISM  list  are  real.  However,  when  dealing  vith  large 
numbers  of  observations,  it  is  always  possible  for  data  to  erroneously  fit  together  to  create  bo¬ 
gus  events.  An  experiment  was  conducted  with  the  help  of  U.SGS  to  see  how  serious  a  problem 
this  might  be.  A  set  of  random  arrival  times  was  generated  for  a  simulated  period  of  ten  days 
with  approximately  the  same  number  of  arrivals  per  day  as  was  reported  during  the  ISM.  These 
data  were  processed  by  the  USGS  program  which  is  normally  used  to  associate  arrivals  and  lo¬ 
cate  events  for  PDE,  Of  course,  all  the  arrivals  in  this  experiment  were  false  whereas  half  or 
more  of  the  arrivals  in  the  actual  ISM  clearly  associate  with  real  events.  Nevertheless,  the 
program  generated  a  surprising  99  bogus  events.  However,  unlike  the  996  ISM  events,  the  99 
events  were  not  subjected  to  careful  analysis,  review,  and  reworking.  For  example,  '7  of  the 
99  events  have  a  RATIO  (defined  above)  larger  than  all  but  two  of  the  actual  ISM  events.  No  ar¬ 
ray  azimuth  and  velocity  information  was  available  for  control.  The  distribution  of  the  stations 
associated  with  an  event  was  not  considered  by  the  program.  Although  it  has  not  been  possible 
to  complete  the  analysis  of  these  99  pseudo-events  without  prejudice,  it  is  our  belief  that  no 
more  than  a  few  would  have  beer  accepted  during  routine  ISM  processing.  Any  accepted  ones 
would  certainly  have  been  graded  either  D  or  DI.  In  summary,  we  conclude  that  no  more  than  a 
small  fraction  of  1  percent  of  ISM  events  might  be  bogus. 

Some  comoarisons  have  been  made  between  ISM  epicenters  and  epicenters  from  other  sources. 
Events  of  grades  D  and  DI  were  excluded  for  these  comparisons.  The  other  sources  of  epicenters 
were  the  U.SGS  PDE  list,  the  LASA  daily  bulletin,  and  the  NORS.AR  daily  bulletin.  The  distance  in 
degrees  from  the  ISM  epicenter  to  each  of  these  was  calculated  and  the  histograms  in  Fig.  I-3(a-c) 
were  generated.  The  histograms  are  constructed  with  boxes  of  unit  width  centered  at  I*.  2*.  . . . ,  20*. 
The  number  of  ISM- PDE,  ISM-LASA,  ISM-NORSAi{  events  used  were  304,  362,  and  285,  respec¬ 
tively.  There  were  5  LASA  and  2  NORSAFl  cases  off  the  scale  of  the  histograms.  These  and 
other  large  discrepancies  are  being  checked  for  specific  explanations  and/or  clerical  errors. 

ISM  and  PDE  epicenters  are  quite  close  to  each  other.  However,  ISM  and  LASA  or  NORSAR  can 
differ  by  several  degrees,  with  NORSAR  showing  more  spread  than  LASA.  This  is  to  be  expected 
since  NORSAR  has  a  smaller  aperture  and  also  gives  bulletin  reports  based  on  PKP  phases.  The 
possibility  that  ISM  epicenters  are  in  error  by  several  degrees  must  be  briefly  entertained,  but 
the  consistency  of  ISM  and  PDE  epicenters  seems  to  make  that  unlikely. 

A  subset  of  32  of  the  almost  200  stations  u«?ed  in  generating  the  ISM  event  list  have  been  se¬ 
lected  for  additional  study.  These  stations  were  selected  to  give  good  global  coverage  while  us¬ 
ing  stations  which  contributed  data  to  large  numbers  of  ISM  events.  Table  1-3  shows  this  list  of 
stations  and  an  indication  of  their  contribution  to  the  ISM  in  terms  of  total  picks  submitted,  total 
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TAHLK  1-2 

DisTHiiu  TioN  OF  ISM  KVFN'rs  :jy  gr\i)i:  and  SFISIMIC  RHGION 


Region 

No. 

Region  Name 

A 

1 

A la.ska- Aleutian  Are 

24 

2 

Fastern  Alaska  to  Vancouver  Island 

1 

5 

California-Novarla  Region 

1 

4 

Raja  California  and  Gu'f  of  t  alifornia 

2 

S 

Mexieo-Guatemala  Area 

9 

6 

Central  America 

4 

7 

Caribbean  Loop 

10 

H 

\ndean  South  America 

17 

9 

ICxtreme  South  America 

1 

10 

Southern  Antilles 

1 

11 

New  Zealand  Region 

12 

Kermadec- Tonga -Samoa  Area 

14 

13 

Fiji  Islands  Area 

5 

14 

New  Hebrides  Islands 

13 

15 

Bismarck  and  Solomon  Islands 

6 

16 

New  Guinea 

4 

17 

Caroline  Islands  to  Guam 

18 

Guam  to  Japan 

57 

19 

Jupan-Kuriles-  Kamchatka 

43 

20 

Southwestern  Japan  and  Ryukyu  Islands 

3 

21 

Taiwan 

4 

22 

Philippines 

8 

23 

Borneo- (’elebes 

1 

24 

Sunda  Arc 

8 

25 

Burma  and  Southeast  Asia 

2 

26 

Ind  ia -Tibet- Szechwan- Yunan 

27 

Southern  Sinkiang  to  Kansu 

28 

Alma-Ata  to  Lake  Baikal 

2 

29 

Western  Asia 

6 

30 

.Middle  Hast- Crimea- Balkans 

12 

31 

Western  Mediterranean  Area 

3 

32 

Atlantic  Ocean 

3 

33 

Indian  Ocean 

1 

34 

Haste rn  North  .America 

t 

35 

Kastern  .South  America 

36 

Northwestern  Europe 

37 

Africa 

38 

Australia 

39 

Pacific  Basin 

40 

Arctic  Zone 

5 

41 

Eastern  Asia 

5 

42 

N.  E.  Asia,  Northern  Alaska  to  Greenland 

1 

43 

Southeastern  and  Antarctic  Pacific 

44 

Galapagos  Area 

1 

45 

Macquarie  Loop 

46 

Andaman  Islands  to  Sumatra 

1 

47 

Baluchistan 

1 

48 

Hindu  Kush  and  Pamir 

7 

49 

Northern  Asia 

50 

Antarctica 

TVBLr  1-3 

CONTRIBUTTOX  OF  32  SELKCTBD  STATIONS  TO  THE  ISM  EVENT  I.IST 
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used  i;.  some  way.  number  of  events  to  which  the  station  contributed,  and  a  breakdown  of  phase 
assignrnP’'.ls.  Of  these  stations,  LAO,  NAO,  HFS,  KBL,  CHG,  and  all  Canadian  stations  were 
^arefjP.y  read  by  an  analyst  or  in  some  other  way  treated  especially  for  the  ISM.  For  tht  other 
stations,  only  the  picks  routinely  suomitted  by  them  to  USGS  for  PDE  were  used  for  ISM.  One 
purpose  i.n  selecting  this  smaller  number  of  stations  was  to  isolate  a  limited  number  of  stations 
for  -nore  detailed  study.  The  other  was  to  configure  a  modest  network  and  evaluate  its  perform¬ 
ance  relative  to  the  overall  ISM. 


TABLE  1-4 

COMPARI.SON  OF  GRADE  DISTRIBUTION  OF  ISM  EVENTS 
WITH  THAT  SUBLIST  GENERATED  AND  GRADED 

USING  ONLY  A  SUBSET  OF  32  STATIONS 

Grade 

Number 
of  ISM  Events 

Number  of  Regraded 

ISM  32  Events 

A 

287 

233 

AI 

61 

67 

B 

56 

0 

BI 

19 

5 

r 

80 

8 

Cl 

32 

8 

D 

197 

239 

DI 

2t4 

268 

An  attempt  has  been  made  to  locate  all  ISM  events  using  only  the  32- station  network.  The 
rt  »nlt  l:?ii  been  encouraging.  It  has  been  possible  to  locate  8i8  of  the  996  events  using  just  the 
32  stations.  Of  the  178  which  could  not  be  lov*ated.  168  of  them  were  in  the  (D,  DI)  category, 

6  in  the  (C.  Cl)  category,  and  4  ‘n  the  B  category  in  the  original  ISM  list.  However,  when  the 
ev»‘nt8  locat-id  with  32  stations  were  graded  there  was  a  modest  downgrading  of  event  quality. 

Tne  number  of  events  in  each  grade  for  the  ISM  and  the  32-station  network  is  shown  in  Table  1-4. 
The  capability  of  smaller  networks  will  continue  to  be  investigated. 

R.  T.  Lacoss 
R.  E.  Needham 
R.  M.  Sheppard 


B.  ANALYSIS  OF  THE  ISM  LONG- PERIOD  DATA 

The  second  half  of  the  ISM  experiment  cDncems  the  detection  and  association  of  surface 
waves  recorded  during  the  month  an  .  the  assignment  of  surface  wave  magnitudes  to  the  events 
located  using  short-peHod  data.  This  will  allow  the  evaluation  of  the  existing  long-period  ar¬ 
rays  and  networks  over  an  extended  period  of  time  using  a  comprehensive  epicenter  list  and 
help  focus  the  problems  that  a  world-wide,  operational,  discrimination  effort  might  encounter. 
The  analysis  of  the  long-period  data  is  not  yet  complete  «»;id  ^\rhat  follows  is  a  brief,  interim 
report. 
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1.  The  Data 


The  long-period  data  for  the  ISM  study  came  from  three  major  sources  and,  for  the  most 
part,  are  based  on  readings  of  seismograms  from  stations  of  various  networks.  The  source 
and  type  of  data  are: 

Department  of  Energy,  Mines,  and  Resources,  Ottawa,  Canada. 

Amplitudes,  periods,  and  associations  of  Rayleigh  waves  recorded 
at  selected  stations  of  the  World-Wide  Standard  Seismographic 
Network  (WWSSN)  and  Canadian  networks  and  the  Hagfors  array, 

Sweden. 

Lamont- Doherty  Geological  Observatory.  Amplitudes,  periods, 
and  approximate  arrival  times  of  Rayleigh  and  Love  waves  re¬ 
corded  at  stations  of  the  Very  Long-Period  Experiment  (VLPE) 
network. 

Lincoln  Laboratory.  Amplitudes,  periods,  associations  and  arrival 
times  of  Rayleigh  waves  recorded  at  selected  stations  of  the  WWSS.N. 

In  addition,  selected  data  from  the  long-period  arrays  (LAS A.  NORSAR,  and  ALPA)  were 
provided  by  the  Seismic  Data  Analysis  Center  (SDAC)  and  Texas  Instrumentr. 

2.  Association 

The  most  difficult  problem  encountered  thus  far  has  been  in  the  development  of  an  optimum 
automatic  scheme  to  be  used  in  the  association  of  the  15,000  long-period  readings  with  the  996 
ISM  events.  Although  in  many  cases  associations  were  made  by  the  readers,  the  criteria  for 
declaring  an  association  were  not  uniformly  set  down  nor  were  the  cases  of  multiple  associations 
adequately  recorded.  The  association  schemes  studied  all  follow  the  same  pattern.  The  velocity 
of  each  Rayleicu  phase  was  computed  using  the  origin  times  of  earthquakes  in  the  previous  two 
hours,  and  this  velocity  was  compared  with  a  theoretical  group  velocity.  If  the  observed  group 
velocities  fell  within  a  window  centered  on  the  theoretical  group  velocity  for  the  appropriate  p>e- 
riod,  an  association  was  declared.  Multiple  associations  were  allowed  and  recorded. 

Initially,  the  theoretical  group  velocity  was  taken  as  an  average  between  pure  continental 
and  oceanic  extremes,  and  velocity  windows  of  ±0.7,  0.5,  and  0.3  km/sec  were  used.  Then  a 
more  elaborate  scheme  to  estimate  the  theoretical  group  velocity  was  employed  in  which  the 
earth  was  divided  into  15*  x  15*  sectors  and  each  sector  was  designated  as  oceanic,  continental.  | 

or  mixed.  For  each  attempted  association,  a  path  was  traced  from  the  station  to  the  hypothetical  J 

epicenter  and  a  path  parameter,  based  on  the  percentage  of  the  path  in  the  various  sectors,  was  J 

I 

assigned.  A  continental  or  oceanic  group  velocity  curve  was  employed  if  the  path  was  more  than  J 

50  percent  either  one.  Otherwise,  an  average  group  velocity  curve  was  used  to  compute  the 
theoretical  group  velocity.  In  this  scheme,  the  observed  velocity  had  to  fall  within  ±0.2  km/scc 
of  the  theoretical  one  in  order  for  an  association  to  be  declare^ 

These  methods  were  tested  on  all  1995  readings  of  main  Ra^*Jigh  phases  made  by  the  Lincoln 
Laboratory  Group  for  the  first  10  days  of  the  ISM.  A  second  test  was  made  on  a  subset  of  those 
readings  with  amplitudes  of  3  mm  or  more  on  the  original  seismogram.  The  results  of  these  . 

tests  are  sho\7n  in  Tables  1-5  and  1-6.  Table  1-5  shows  the  number  of  times  that  all  individual  i 

readings  and  those  of  the  subset  were  associated  with  an  event  using  the  various  association 
methods.  In  the  sense  that  the  method  which  yields  the  maximum  number  of  unique  associations 
is  the  optimum  one,  an  average  group  velocity  curve  with  a  0.5- km/sec  window  would  appear  to 
yield  the  best  results  of  all  the  methods  tested  when  all  tlie  readings  are  considered.  The  per¬ 
centage  of  unique  associations  increases  for  all  methods  when  only  readings  of  amplitude  >  3.0  mm  I, 
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TABLE  1-5 


NUMBER  OF  TIMES  THAT  LONG-PERIOD  READINGS  MADE 
DURING  FIPST  iO  DAYS  OF  ISM  WERE  ASSOC’I  \TED  WITH  EVENTS 


Number/percentage 
of  all  readings  (19Q5) 

Number  of  Associations 
of  a  Single  Reading 

Method  of  Association 

0 

1 

2 

3 

a 

793 

39 

537 

32 

968 

48 

1184 

59 

878 

44 

9  36 
47 

841 

42 

741 

38 

1 

H 

8 

0.4 

11 

0.6 

4 

0.2 

2 

0.1 

Average  velocity 
±0.7  km /sec  window 

Average  velocity 
±0.5  km /sec  window 

Average  velocity 
±0.3  km/sec  window 

Path-dependent  velocity 
±0.2  km /sec  window 

66 

184 

90 

23 

7 

Average  velocity 

17 

49 

24 

6 

2 

±0.7  km/sec  window 

92 

199 

67 

12 

2 

Average  velocity 

Numher/percentage 

24 

53 

18 

5 

0.5 

±0.5  km /sec  window 

of  readings  ^3.0  mm 

(372) 

149 

182 

38 

2 

1 

Average  velocity 

40 

49 

10 

0.5 

0.3 

±0.3  km/sec  window 

156 

196 

19 

1 

0 

Path-dependent  velocity 

L 

42 

53 

5 

0.3 

0 

±0.2  km /sec  window 
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TARI.L  1-6 

MMni.'H  AND  LKfiCIINTAC.L.S  OF  AG  HLFAIFNTS 

ri:t\vi:i:n  altomxtk'  association  sciilmks  and  trosl 

M\Di:  HV  IILND  WHLN  ONLY  ONK  AFTOMATK' 
ASSOCIATION  WAS  .AlADi: 

All 

Like  Associations 

Method 

Total  Single  Associations 

Percentage 

Average  velocity 

604 

0.7° 

±0.7  kni  /sec  window 

878 

Average  velocity 

740 

0.79 

±0.5  km /sec  window 

936 

.Average  velocity 

664 

0.79 

±  0. 3  km/sec  window 

841 

Path-dependent 

61  3 

0.82 

velocity  ±  0.2  km /sec 

741 

window 

.Amplitude  3.0  mm 

Average  velocity 

155 

0.84 

±0.7  kni/sec  window 

184 

Average  velocity 

164 

0.82 

±0.5  iun /sec  window 

190 

Average  velocity 

147 

0.80 

±  0. 3  km /sec  vvindow 

i  82 

Path-dependent 

170 

0.86 

velocity  ±0.2  km /sec 

196 

window 

are  used.  This  increase  is  most  <lramatic  for  the  method  which  used  a  path-dependent  group 
velocity  and  a  smaller  association  window.  Assuming  that  the  period  and  arrival  time  of  the 
larger  phases  were  read  with  greater  accuracy  than  smaller  ones,  we  must  conclude  that  path- 
dependent  velocities  are  justified  when  the  data  were  based  on  large,  clear  phases. 

Table  1-6  summarizes  another  test  that  was  made  of  the  association  algorithms  using  the 
associations  made  by  the  readers  at  Lincoln  Laboratory.  Rased  on  the  same  readings  as 
Table  1-5,  this  table  shows  the  number  and  percentages  of  times  that  a  unique  automatic  associ¬ 
ation  agreed  with  the  Lincoln  readers'  association.  If  the  latter  are  assumed  to  be  the  correct 
association,  then  the  various  automatic  methods  can  be  tested  on  the  criteria  of  number  of  cor¬ 
rect  unique  associations  and  percentage  of  correct  associations  based  on  the  number  of  unique 
associations.  Table  1-6  shows  that  w'hen  all  the  data  are  considered,  the  average  group  velocity 
curve  with  a  0.5-km/sec  window  gives  the  maximum  number  of  agreements,  although  the  per¬ 
centage  of  like  associations  made  by  this  method  is  the  same  or  less  than  the  others.  When  the 
larger  amplitude  data  are  considered,  the  method  which  uses  the  path-dep.?ndent  velocities  is 
slightly  better  than  the  others  tested.  Clearly,  criteria  other  than  group  velocity  could  be  used 
in  the  automatic  association  of  surface  wave  readings  with  a  list  of  earthquakes.  The  amplitude 
of  the  surface  wave  considered  with  the  body  wave  magnitude,  depth,  and  distance  of  the  event 
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would  certainly  help  sort  out  multiple  associations.  However,  such  a  scheme  would  have  to  be 

applied  with  care  in  a  discrimination  experiment. 

Although  work  continues  on  an  optimum  surface  wave  association  algorithm,  we  have  applied 

the  one  based  on  an  average  group  velocity  curve  and  a  ±0. 5-km/sec  window  to  the  entire  ISM 

long-period  data  set.  Using  only  readings  which  possess  a  unique  a.ssociat  o.-'.  and  the  surface 

v/ave  magnitude  formula  M  =  log(A/T)  +  1.66  log  A  +  3. 3  (a  vertic?.!  Rayleigh  amplitude  in  mi- 

s 

crons,  20“  <  A<  160“,  and  15  <  T<  30  sec),  we  have  computed  an  average  M^  based  on  two  or 
more  readings  for  264  of  the  827  ISM  events  to  which  a  body  wave  magnitude  (mj^)  was  assigned. 

In  Fig.  1-4,  we  inevitably  plot  this  average  M  vs  the  average  rr..  of  these  events.  Although  there 
are  yet  errors  in  association  represented  in  Fig.  1-4,  it  gives  us  something  to  test  subsequent 
schemes  against.  In  order  to  be  considered  improvements,  jther  methods  of  magnitude  deter¬ 
mination  or  reading  association  should  reduce  the  scatter  ir.  or  add  more  points  to  Fig.  1-4,  or 
both.  Figure  1-4  does  produce  one  encouraging  note;  the  circled  point  is  a  presumed  underground 
nuclear  test  at  the  Semipalatinsk  test  site  in  Eastern  Kazakh.  Without  qualifying  the  remaining 
events  by  depth  or  region,  this  presumed  explosion  separates  from  them  rather  well. 

J.  Filson 
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Fig.  I-l.  Incremental  histogram  of  ISM  m^j  (AVEMB)  for  (a)  events  graded, 
A,  AI,  B,  BI,  C,  and  Cl,  and  (b)  events  graded  D  or  DI  which  have  an  as¬ 
signed  magnitude. 


Fig.  1-2,  Geographic  distribution  of  ISM  events-  (a)  all  events,  and  (b)  events 
of  grades  A.  AI.  B.  BI.  C.  and  Cl. 
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NUMBER  or  EVENTS 


Fig.  1-5.  Ili.stograms  comparing  A,  AI,  f5,  151,  C,  and  Cl  grade  ISM  epicenters 
with  (a)  r.SCIS  Preliminary  Determination  of  Epicenters  list,  (b)  I..ASA  dail\ 
bulletin,  and  (c)  N()R.SAR  dailv  bulletin.  I  nits  are  degrees. 


AVERAGE  m. 


Fig.  1-4,  Average  Mg  vs  average  for  264  events  of  ISM.  Average 
determinations  are  based  on  two  or  more  measurements  in  restricted  pe¬ 
riod  and  distance  range.  Only  unique  associations  from  one  of  association 
schemes  was  used.  Circled  point  is  presumed  explosion  in  Eastern  Kazakh. 
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II.  AHIIAV  S1‘l  l)!i:S 


A.  NOKSAK  LIPlCENTKIi  AZIMI  TH  KI^IIOIJS 

A  continuing  effort  to  understand  the  XOKSAK  «'picentt*r  anomalies  is  producing  more 
intare-.ting  results  than  tlie  I. ASA  epicenter  anomalies.  The  NOHSAlt  azimuth  anomaly,  in  par¬ 
ticular,  is  one  that  is  difficult  to  understand  or  approximate  with  a  model.  'I'he  azimuth  anomaly 
is  defined  as  the  difference  between  the  true  azimuth  to  tlu  epicentei-  and  the  actual  olisetwed 
azimuth  at  NOItSAK,  A  least-square  fit  was  made  to  the  subarray  arrival  times  obtained  from 
the  event  processor  that  were  corr'ected  to  remove  any  misalignment.  I'rom  this  process  the 
velocity,  azimuth,  and  epieent<‘r  were  obtained.  When  the  arrav  epii cntei'  is  plotted  with  the 
true  PDK  epicenter  (Fig.  II- 1),  the  azimuth  error  is  not  finlv  appnia-nt  hut  is  <]Uite  profound.  Be¬ 
cause  of  the  character  of  this  anomaly,  the  primary  cause  must  he  the  array,  with  minor  sec¬ 
ondary  effects  caused  by  the  epicenter.^  One  cannot  account  for  the  azimuth  error  as  a  result 
of  a  plane-dipping  interfaee  at  the  MOIK).  Only  in  a  ver  v  general  way  would  a  plane-dipping 
Mono  approximate  the  azimuth  or  dT/dA  error.  .More  likely,  one  needs  a  dipping  MOIIO  with 
an  egg-shaped  bulge  under  NOHSAH.  The  dip  of  the  .MOIIO  would  be  westerly  while  the  axis  of 
the  dome  in  the  MOHPJ  would  be  west-northwesterly.  Kay  tracing  is  now  being  used  to  model 
the  actual  data. 

H. M.  Sheppard 

B.  ANALYSIS  OF  AKKAV  DlAOKA.MS  A  T  LASA  AND  NOKSAK 

A  powerful  method  for  the  analvsis  of  approai  h  angle,  or  slowness,  anomalies  for  P-wave 
signals  at  an  arrav  is  ha.sed  on  their  representation  in  .s|own<*ss  hlT/dA,  azimuth)  spare.  Such 
a  representation,  known  as  an  array  niagram,  involves  the  nlotting  of  these  anomalic.s  a.s  arrows, 
in  which  the  tail  represents  the  dT/dA  and  azimuth  measurt'd  for  the  P-vvave  signal  at  the  array, 
and  the  head  represents  the  dT/dA  and  azimuth  that  would  be  obtained  from  the  true,  or  in  this 

case  NOS,  location  and  J-B  tables.  This  approach  has  been  used  by  Manchec  and  W’eichert,^ 
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Vinnik  and  Nikolayev,  and  by  Davie.s  and  Slu'ppard,  and  C  apon  at  LASA. 

A  problem  of  considerable  importance  is  tbe  determination  of  the  various  inhomogeneous 
regions  within  the  earth  that  are  responsible  tor  these  f^-wave  slowness  anomalies  measured 
by  LASA.  These  inhomogeneities  may  he  described  broadly  as  due  to  the  following  causes: 

<1)  Kandom  part  of  local  crust  and  upper- mantle  structure  at  the  receiver, 

(2;  Deterministic  part  of  local  crustal  structure  at  the  receiver, 

(3)  Propagation  path  effects,  primarily  in  the  lower  mantle,  and 

(4)  Source  region  effects,  due  to  structures  such  as  dipping  plates  of  lithosphere. 

It  has  be°n  shown  by  (!apon^  that  the  eifect,  due  to  the  structure  listed  in  item  (1 ),  on  the 
P-wave  slowness  anomaly  is  negligible.  The  reason  for  this  is  that  the  LASA  tends  to  average 
out  the  effecl  due  to  this  .‘■tructurc  when  the  slowness  of  the  P-wave  is  measured  by  means  of  a 
least- squares  fit  of  a  plane  wave  set  of  arrival  times  to  the  actual  set  of  arrival  times  measured 
over  the  entire  LASA.  This  means  that  the  source  of  the  anomaly  must  then  be  due  to  some,  if 
not  all,  of  the  structures  listed  in  items  (2),  (3),  and  (4). 

If  the  anomaly  were  caused  by  a  plane-dipping  interface  directly  under  the  array,  the  array 
diagram  would  consist  of  a  set  of  error  vectors  whose  magnitudes  were  all  approximately  the 


s.trnc  ;in(l  nil  pointfd  prcf  isrl\  in  tliv  (i i tn-i  t dp  of  dip  of  ilif  interface.  Thus,  the  effect 

of  a  sl'ipiivu  iiitorf.'uc  can  he  ol)ia!iu'd  from  tlv  arra\  diauram  hv  littine  a  ia>nsinnt  vector  to  the 
error-  \eetoi-s  in  the  diacram.  These  i-esiihs  •.vei'e  first  shown  hy  .Manchee  and  Weichert,  and 
an  example  of  this  has  been  uiveii  bv  (  apon  Dav  ies  and  .Sheppard  havi*  pointed  out  that  a  slowly 
varvinti  (list ribiit ion  of  arrows,  all  approximatidv  of  euLUil  lon^th  and  pointi'*^*  in  the  same  direc¬ 
tion,  would  indicate  a  bias  mtroducecl  beneath  the  array.  This  bitis  may  not  be  due  necessarily 
to  a  dipptne  inti-rfnce  but,  for  example,  from  ;i  distribution  or  coating  of  incremental  time  de¬ 
lays  that  are  not  independent  of  eaeb  other  foe  the  various  directions  of  approach.  That  is,  these 
incremental  time  delavs  are  di reel  i  >n;ill>  independent.  Thus,  the  effect  of  the  deterministic 
part  of  tile  local  criistal  striicture  at  the  array,  tii\(Mi  in  itmn  (2), may  be  obtained  by  detertnining 
the  averaye  of  all  the  error  vectors  in  the  arr.av  diagram.  l  i>r  simplicitv,  this  structure  may 
be  e()nsidered  to  be  represented  by  some  eipiivaletit  plane- dipping  interface,  and  the  dip  and 
strike  direction  for  it  may  be  obtained  easily  from  the  mat'iittiide  and  direction  of  the  average 
vector  bv  using  some  results  due  to  Nia/t.’  Once  the  effet  t  due  to  this  .Structure  has  been  re¬ 
moved,  by  subtracting  the  average  of  all  the  erroi  vectors  from  on  h  of  the  vectors  in  the  array 
diagram,  then  the  source  of  the  remaining  anomaly  must  be  due  to  the  structures  in  items  13) 
and  (4). 

In  the  previous  .SA’l'S,^  a  method  of  nntil.v.sis  based  on  the  array  diagram  was  presented  which 
showed  that  the  crust,  upper,  and  mid- mantle  structures  directly  beneath  LASA  were  responsible 
for  most  of  the  P-vvave  slowness  aiionialv.  The  proof  was,  of  necessity,  only  heuristic  in  nature 
since  there  is  no  unujiie  wav  m  which  earth  structures  can  he  defined  as  ritsponsible  for  the  anom¬ 
aly.  However,  additional  results  are  now  presented  to  support  the  proof  presented  previously, 
and  similar  rt’sults  for  the  .NOUSAK  are  also  presented. 

The  basis  of  the  proof  prt'Stmted  in  the  prev  ious  .S A'l'.S  ^  is  that  the  .spatial  rate  of  variatifin 
of  the  error  vectors  can.  at  least  form.ullv,  be  asso.-inied  with  heterogeneities  at  certain  depths 
in  the  earth.  The  more  rapid  thi.se  sp:iiial  v.-iriations  hecoim',  the  more  distant  from  L.ASA  thi 
heterogeneities  must  be  placed  to  explain  the  .'inomal ie.s .  Some  artificial  array  diagrams  are 
now  presented,  for  hypothetical  earth  structures,  to  support  this  eoi  dusion,  .All  these  results 
were  obtained  by  using  a  ray- tracing  prog’  am  vvliieh  is  based  on  the  computational  procedure  of 
Nia/.i^^  for  finding  the  slowness  of  a  rjiy  which  is  refrtieted  hv  a  plane-dipping  interface.  The 
first  example  deals  with  the  corrugated  surface  proposed  by  Croenfield  and  .Sheppard,^  as  a  model 
for  the  Mt)ll()  under  LAS.A,  in  order  to  explain  the  source  of  the  anomaly.  The  array  diagram 
for  this  structure  is  shown  in  I'lg.  11-2,  and  it  is  seen  that  a  relatively  slowly  varying  distribution 
of  arrows,  all  approximately  of  eiiiial  length  and  pointing  in  almost  the  same  direction,  is  ob¬ 
tained.  This  result  agrees  with  the  notion  that  such  a  distribution  of  arrows  should  he  associated 
with  a  struv  turc  located  l  enenth  the  array. 

The  average  vector  in  F’ig.  11-2  has  a  magnitude  of  O.ob  sec/deg,  and  points  at  an  azimuth 
of  34C“.  By  using  the  pr  u  edure  of  Niazi^  once  again,  it  is  possible  to  show  that  this  error  vec¬ 
tor  is  due  to  a  hypothetical  plane-dipping  interface  under  LASA  with  dip  of  about  H*  and  direction 
of  dip  at  an  azimuth  of  340”,  assuming  that  the  ratio  of  crustal-to- mantle  velocity  is  0,75.  These 
parameters  corresp(jnd  roughly  with  those  of  a  plane  which  might  be  considered  as  a  least- squares 
fit  for  the  corrugated  surface.  That  is,  this  average  plane  can  be  considered  as  one  which  inter¬ 
sects  the  corrugated  surface  in  such  a  way  as  to  minimize  the  sum  of  the  squares  of  the  deviations 
between  them.  The  error  vectors  in  the  array  diagram  in  Kig.  11-2  can  be  considered  as  due 
primarily  t(^  this  average  plane,  with  small  deviations  introduced  due  to  the  effect  of  the  corru¬ 
gated  surface. 


16 


'riu  riccond  oxamph-  deals  u  itn  a  pyt-amidal  struetufe  ui'ose  apex  is  located  I'iO  km  beneath 
the  (-(Miter  of  LASA.  'flie  dip  an^fle,  for  all  four  sides  of  the  pyramid,  is  10°  and  the  dip  direction 
of  OIK*  of  the  sid(’S  is  at  an  azimuth  of  45",  so  that,  of  eoursc,  the  clip  directions  for  the  other 
sides  are  i  i=>°.  Z25",  and  31  S".  The  velocity  contrast  across  the  striieture  taken  as  O.SS, 

The  array  diagram  for  this  structure  is  shown  in  h’iy.  11-3,  and  it  is  seen  that  a  relativelv  rapidly 
varvin^  distribution  of  arrows  is  obtained.  It  is  recalled  from  l-’i^e  111-6  of  th>;  previous  SAT 
that  the  horiztJiital  cross  section  of  the  disk  in  the  75-  to  225-km  depth  range,  at  a  cleptli  of  about 
150  km,  is  e(jual  to  about  four  tiiiies  the  area  of  the  LASA  aperture.  As  a  consequence,  it  was 
shown  that  there  should  he  four  independent  err(»r  vectors  introduced  in  the  array  diagram  by 
a  struetuiM'  located  at  a  depth  of  I kni.  it  is  expected,  therefore,  that  the  array  diagram  in 
l-'ig.  ll-  }  should  contain  four  indc'per.dent  error  vectors.  This  is  not  exactl/  true,  but  can  he 
considered  as  a  reasonable  apjo-oxi mation.  !-'oi*  example,  the  array  diagram  in  Fig.  11-3  can  he 
divided  into  four  equal  (piadrants,  and  within  i-ach  particular  quadrant  the  average  ••'ictor  can 
he  found,  and  it  can  he  ust_d  to  riqilace  the  vectors  within  the  quadrant.  The  array  diagram  ob¬ 
tained  m  this  manner  contains  four  independent  vectors  and  can  be  considered  a  good  apnrt'xi- 
mation  for  the  re.sult  in  Fig.  11-  3.  it  should  he  noted  that  the  orientation  of  the  axes  used  to  ob¬ 
tain  the  (jundrants  is  no  jiarticulnrly  important,  as  a  g(  od  approximation  is  ohtaint'd  for  all  such 
ori(*ntations.  Thus,  this  examph'  also  provides  further  corroboration  for  the  heuris  ic  proof 
pre.sented  in  the  last  S.A'I'S. 

Ihe  analysis  procedure  described  previouslv^  has  been  applied  to  the  array  diagram  at 
.'iORSAIi.  't  he  decomposition  of  the  region  direetly  under  .NOKSAK  is  shown  in  Fig.  11-4,  and 
differs  .somewhat  from  the  corresp'Muling  result  for  LASA  shov\  n  in  I-'ig.IlI-6  of  the  pr  evious 
.SA’IS,  since  th<'  aperture  of  \()H.SAl{  is  110  km  as  compared  with  200  km  for  LASA.  The  orig- 
iiuil  data  consisttMl  of  7  22  (’rror  vec'tors.  However,  in  order  to  provide  uniform  Cviverage  in 
the  array  diagram,  a  data  decimation  s'-heme  was  used  in  which  the  array  diagram  was  divided 
into  a  grid  of  20  <  20  cells,  where  each  cell  is  a  square  of  1  x  1  sec/deg.  The  first  ten  vectors 
to  fall  within  a  giv«Mi  .  ell  were  taken  as  acceptable  data.  'I  his  led  to  a  total  of  218  error  vectors 
as  shown  m  Fig.  11-5.  The  a\erage  of  these  error  vectors  was  computed  and  found  to  point  in 
the  direction  corresponding  to  an  azimuth  of  2H6°,  with  a  modulus  of  0.52  sec/deg.  This  mean 
vector  could  he  associated  with  a  plane-dipping  interface  which  dips  down  at  an  azimuth  of  286°. 
These  pnrami’ters  could  he  associati’d  with  a  mean  dip  of  6°  of  the  .MOHO  under  NOKSAR,  as¬ 
suming  a  crustal- to- mantle  velocity  ratio  of  0.75,  although  the  MOHO  is  not  necessarily  the 
only  source  of  anomaly. 

In  l-'ig.  11-6  we  have  tin-  sum  of  the  average  vector  and  all  the  vectors  obtained  using  the 
suh-a veraging  method  descriheti  pn-vioiislv.^  The  great  similarity  of  the  array  diagram  in 
Fig.  11-6  with  thac  in  l-'ig.  11-5  should  he  noted.  This  similarity  is  also  borne  out  by  the  small 
magnitudes  of  the  v(‘ctoi-.s  shown  in  Fig.  11-7  which  depicts  the  difference  between  the  error  vec¬ 
tors  in  l-‘igs.  11-5  and  II-6.  The  mean-sr|uare  value  of  the  error  vectors  shown  in  Fig.  il-7.  rel¬ 
ative  to  that  in  the  original  arrav  diagr."m  shown  in  !-’ig.  11-5,  is  about -13.5  dB.  This  means 
•  i.at  the  analvsis  procedure  has  been  able  to  remove  more  than  96  percent  of  the  original  error 
vectors,  when  measured  in  terms  of  mean-sejuare  value  of  the  magnitudes.  This  should  be 
compared  with  the  corresponding  result  for  I.ASA  of  90  percent  obtained  previously. 

It  is  also  possible  to  find  the  mean-squa^’e  value  of  the  error  vectors  in  each  depth  range, 
relative  to  that  in  *he  original  arras  diagram  shown  in  Fig.  II- 5.  The  rcsult.s  of  such  a  compu¬ 
tation  are  shown  in  I-‘ig.  11-8,  where  the  ordinate  has  been  labeled  as  relative  amount  of  hetero¬ 
geneity.  It  is  seen  from  this  figure  that  there  is  an  approximately  decreasing  amount  of 
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heterogeneity  as  the  depth  inc.reases.  This  agrees  roughly  with  tlic  notion  that  the  eai  th  l)(.romes 
more  homogeneous  as  the  depth  increases,  and  agrees  reasonably  well  with  the  ('orrcsponding 
result  obtained  previously  for  LASA.^ 

It  should  be  mentioned  that  the  orientation  of  the  axes  used  te  divide  the  arrav  diagram  into 
cells,  for  the  purpose  of  obtaining  the  sub-averages,  is  allowed  to  vary  in  six  steps,  from  0° 
to  7  5“,  with  increments  of  i5®  between  the  steps.  The  results  for  the  various  orientatuins  are 
averaged  to  obtain  th^  Tnal  sub-averages.  In  addition,  ttie  results  for  the  various  orientation 
angles  ere  similai  co  each  other.  This  should  sei’ve  to  eliminate  any  criticism  the  analysis 
procedure  is  sensitive  to  the  orientation  angles  employed  to  obtain  the  sub-averages.  Another 
refinement  used  in  the  method  is  to  apply  a  two-dimensional  triangular  weighting  funriion  to  the 
data  in  each  cell  when  the  sub-averages  are  obtained.  This  is  done  to  prevent  the  analysis  method 
from  being  sensitive  to  large  error  vectors  near  the  boundaries  of  a  cell.  However,  the  results 
obtained  either  with,  or  without,  this  type  of  weighting  are  quite  similar. 

It  is  also  possible  to  sum  the  average  vector  and  all  the  vectors  obtained  using  the  sub- 
averaging  method  at  points  in  the  array  diagram  which  do  not  necessarily  represent  the  original 
data  points,  e.g.,  on  a  uniform  grid.  This  ha.*!  been  done  at  both  NORS.AR  and  L.ASA  and  the  re¬ 
sults  are  shown  in  Figs.  II- 9  and  11- 10,  respectively.  These  data  can  serve  to  define  the  P-wave 
anomalies,  and  thus  the  station  corrections,  from  aseismic  regions.  It  should  be  noted  that  the 
anomalies  to  be  expected  for  the  core  phases  at  NORSAR  and  L.AS.A  are  also  depicted  in  Figs.  11-9 
and  II- U  respectively. 

.1.  Capon 

C.  A  RANL»OM  MEDIUM  ANALYSIS  OF  CRUST  AND  UPPER- MANTLE 

STRUCTURE  UNDER  NORSAR 

The  measurement  of  the  P-wave  amplitude,  phase,  and  slowness  within  subarrays  at  LASA 

has  been  considered  extensively  in  the  past.  These  analyses  have  revealed  substantial  variations 

in  these  quantities  between  subarrays.  Until  recently,  there  was  no  adequate  gcophvsical  mode! 

of  the  crust  and  upper  mantle  under  LASA  which  could  explain  the.se  P-wav€  anomalies.  However, 
8  5 

recent  work  by  Aki  and  Capon  has  led  to  the  development  of  earth  models  which  could  success¬ 
fully  explain  the  source  of  these  anomalies  in  certain  cases.  This  recent  work  is  based  on  the 

notion  of  characterizing  the  crust  and  upper- mantle  structure  under  LASA  as  a  random,  or 
9 

Chernov- type,  medium. 

A  similar  analysis  has  now  been  done  at  NORSAR.  Tow'ard  this  end,  the  P-wave  amplitude 
and  phap'»  spectra  have  been  measured  for  a  total  of  32  events  distributed  at  various  azimuth.s 
and  distances  from  NORSAR.  These  spectra  were  measured  over  a  5-sec  interval  starting  at  a 
time  which  is  determined  by  the  arrival  time  at  a  sensor  of  the  best-fitting  plane  wave  propagating 
across  NORSAR.  The  slowness  of  this  plane  wave  is  determined  from  the  P-wave  beamforming 

g 

location  system  at  NORSAR.  This  procedure  is  similar  to  that  used  by  Aki.  The  body-wave 
magnitudes  of  the  events  employed  in  the  study  were  all  larger  than  about  5.0,  in  order  to  insure 
a  reasonably  good  s:gnal-to-noise  ratio. 

There  were  9  subarrays  and  6  sensors  per  subarray  at  which  spectral  m.ea.«iurc-ments  were 
made.  'Fhese  54  measurements  were  divided  into  three  sets,  and  within  each  set  13  sensors 
were  considered.  The  subarrt  /s  used  in  each  set  are: 

Set  Subarrays 

1  Oin,  02B,  03H 

2  OlA,  043,  06C 

3  0513,  06B,  07 R 
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These  sets  of  subarrays  were  chosen  so  that  wii  lin  each  set  the  subarrays  are  located  on  similar 

geologic  struct  "  'S  at  NORSAR.  The  particula’  frequencies  used  in  the  analysis  were  0.6,  0.8, 

1.2,  1.6,  and  2.0  Hz.  These  five  frequencier  were  chosen  since  they  are  similar  to  those  used 

by  Aki^  and  Capon"^  in  their  related  work  at  LASA. 

8  5 

The  method  which  has  been  used  by  Aki  and  Capon  to  test  whether  a  medium  can  be  con- 

o 

sidered  to  be  randomly  inhomogeneous,  in  the  sense  given  by  Chernov,  is  as  follows.  The 

9  2  2 

theoretical  relationship,  due  to  Chernov,  between  B  /S  and  R^^  for  various  values  of  the  wave 

parameter  is  plotted,  as  shown  in  P’ig.  P**!!,  where  and  are  the  mean-square  fluctuations 
of  amplitude  and  phase,  respectively,  and  is  the  crosscorrelation  coefficient  between  the 
amplitude  and  phase  fluctuations.  The  90- percent  confidence  region  is  also  shown  in  Fig.  11-11, 
as  determined  from  the  size  of  the  sample  used  i  i  the  investigation  The  method  is  to  estimate 

B^S^  and  R^^  from  the  P-wave  amplitude  and  phase  fluctuations  and  determine  whether  these 
estimates  fall  within  the  90-percent  cc  nf*dcnce  region  predicted  by  the  Chernov^  theory,  as  given 
in  Fig.  11-11.  These  estimates  have  been  made  at  NORSAR  for  the  five  frequencies  and  the  three 
sets  of  sensors  indicated  previously,  and  the  results  of  these  15  measurements  are  shown  in 
Fig.  11-11.  It  is  quite  clear  from  the  figure  that  all,  except  one,  of  the  measurements  fall  out¬ 
side  the  confidence  limits.  This  indicates  that  the  crust  and  upper- mantle  structure  under 
NORSAR  does  not  appear  to  satisfy  the  conditions  required  for  a  Chernov-type  medium. 

Measurements  were  also  made  of  the  transverse  autocorrelation  coefficients  of  the  ampli¬ 
tude  and  phase  fluctuations  at  78  (=13*  12/2)  spatial  lags  and  six  successive  values  were  averaged 
together,  to  improve  the  stability  of  the  results,  sc  that  a  total  of  13  points  was  obtained.  These 
measurements  are  shown  in  Figs.  11-12  and  11-13  for  a  frequency  of  1.2  Hz  and  for  the  second 

Q 

set  of  senbors,  along  with  the  theorct’eal  predictions  clue  to  Chernov.  These  results  are  typical 
of  those  obtained  at  other  frequencies  and  fer  other  sets  of  sensors.  The  results  i''  Figs.  11-12 
and  11-13  show  quite  clearly  that  the  correlation  distance  for  this  particular  set  of  measurements 
is  about  4  km.  However,  if  the  entire  set  of  15  measurements  is  considered,  it  is  found  that 
the  correlation  -lis^ance  varies  between  about  4  and  8  km,  with  many  of  the  measurements  be¬ 
tween  4  and  6  km.  Thus,  it  appears  that  the  average  correlation  distance  for  the  15  measure¬ 
ments  is  about  6  km. 

The  interpretation  of  these  results  is  as  follows.  The  wavelength  of  the  P-wave  in  the  0.6- 
to  2.0-Hz  range  varies  from  about  13  to  4  km,  assuming  an  average  compressional  velocity  of 
8  km/sec  in  the  upp>er  mant)e.  Thus,  the  wavelengths  are  larger  or,  at  best,  of  the  same  order 
as  the  size  of  the  inhomogeneity,  or  correlation  distance.  This  violates  the  condition  of  large- 
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scale  inhomogeneities,  assumed  by  Chernov,  whi-:h  requires  that  the  wavelength  be  much  smaller 
than  the  correlation  distance,  lienee,  it  appears  that  the  Chernov  theory  may  not  be  applicable 
to  the  probleir.  of  the  scatterii.g  of  the  P-waves  under  NHRSAR.  In  particular,  it  appears  that 
multiple,  and  not  single,  scatter  i nr  is  taking  place  under  NORSAR,  so  that  the  conditions  re¬ 
quired  in  the  soluifon  of  the  wav  «"^uation  by  means  of  the  Born  approximation  are  not  satisfied. 

0 

This  conclusion  was  also  reached  by  Aki  in  his  analysis  at  LASA  for  frequencies  greater  than 
1  Hz. 

J.  Capwn 

K. A.  Bertcussen 
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1).  AMl’LlTl  l)i:  AM)  I’llASi;  \  AIMA  1  H  >N.s  ( )I  MiAAl.s 

WITHIN  NOHSAH  AND  I. ASA  Si  HAilHAiS 

P-\vnvf  nmpliludr  and  |)Ii;«s«  ihu-iuati  'iis  nif  ijt-UiL  >?ii'i;fd  ,i!  ^  !•■,>,  .  -  \  end 

LASA.  'Ihc  an’i{)l  i  tudc  ri  iH  t  iin !  I' tiii-  m  isni  'Ui  n  ni.s  if*  iMicej  n  'i|i|  i  r  >  I  •  m  ir* 

snvtTf  tl'an  those  at  I, ASA,  the  inam  i  t*  ■>  >n  i-fitij  Hi  n  Hn'  In  -Ii  'i  .  hi  ^ifin  MUf  >111.-:.  r  \<MA\1 
tran'sniits  a  wide  hand  oT  i t  t  i|iieni’ les  <  si»fi  lalb.  I'n  ate  1  ti,nn  1  11/  tin  m  ii-  n  ■;  11  iv.ev.  i  . 

LASA  IS  located  on  crustal  sftlitin-tiis  ovee  a  '  iv\-L’  ni'iii  i  ntanih  .  v.liiefi  atti  irrde.'.  hu.'h 
fi’cquencies. 

As  a  function  of  fre(|iK-ncv,  du  anii)lit'id«  'Ln  Hi  it  nii';  u  itietii  snhro-i  e,  S  ti  e  sui- 

prisin^lv  similar  to  ilio>e  witlim  I.ASA  sul-an  avs  II>  i  *  .■  >  pr-  .1  n'  :■  'i.e  >>list  r.ati  in  'f  nm- 

plitudc  and  phase  variation  wnhm  tin-  sni  la  rra  .  •:  NolAAl.  .tin.  !  's  \  In-  ,1  pi  •  -‘nued  .■pi  .ion 

in  Eastci’n  Kazakh  on  h  .lune  !  .»7 1 .  in  oile  r  t  >  eipnli/t  Hi.  il  .i  r.i  s’..  .n  •  .•  '  n  1  :  .  >1  . 

outer  rin^  sensors  and  tht*  centei’  sensor  of  e.i.  li  snltai  r;!\  at  I.A.sN  v.ei-.  u>.  d.  reid  -i  NidpsM; 

all  five  outer  subarray  sensors  and  center  sensors  ut  re  u  d  In!  oa  I!-  I  )  .-h  >u,;  1  .nip, ins  ui 

between  typical  subarrays  at  LASA  and  NOHSAIl. 

At  each  subarray,  the  l-oiiriet  speeii-al  ratio  of  eael;  m-hs  n*  to  (}(•  sul  an  as  be  ne  v..,;-; 
I'omputcd  usin^  10  sec  tjf  P-wave  data,  from  the  spemral  ratio;;,  ?!,■  hk  m  r.nd  staiui.  nl  flesi 
a'lSon  of  the  amplitude  and  phase  didav  spi'ctro  sseia*  e  imputed  f  n  •  at  h  snh-.i  r:»s.  1  !u  pii.ise 

delay  i.s  defined  here  to  be  the  time  del.av  in  seconds  •!  eat  h  fret  •  ae',  t*  nop  .le  m.  i*  c  l,i'v‘  n 
by  -<(>{a‘)/u;,  where  u?  is  in  radians /second  and  ,,(u’ 1  is  tit-  nnvrapped  phase  "f  the  ^p  (  ir.il  r.ifiti. 

Figure  H-l*'  shows  the  mean  atui  standard  deviation  u  tin  aUiplinidi  sfc  <  tral  rrPi..  it  tvo 
subarrays  at  NOKS.Mi  and  LASA.  Subarravs  «>(’  at  NOlt.SAH  .ind  I A  m  I.AsA  ire  ivpanl  oi  ito- 
best  quality  subarrays  for  this  event.  Tlu'  nu  an  .nmplituile  sf.eetral  rniio  of  si  nste  to  siil.aria\ 

beam  is  very  flat  over  the  freiiueney  band  I  to  2  A  H/.  tor  both  suhnrravs.  Ilu  da  sin  1  hiu  s  siiow 

die  amplitude  ratio  plus  and  minus  the  standard  deviation  .let  ,ss  t  .i.  !i  sid>arra\  1  h  |  .u  Hm  ;.* 
plots  are  the  mean  phasi’  tlelnv  in  seconds  with  standard  ;!•  viaii  'H  pter  Imnts  foi-  tip  srnn* 
subarravs.  l-’or  bt»th  subarravs,  the  mean  plnase  dei.'.v  is  n.'arlv  n  .'-i-.  indie  1' mg  ih.it  .ai!  fi  i 
qiiencv  components  of  the  single  stmsoc  an  t  1  pliase  witli  lUr  beain  eompoueiu.;. 

In  Fig.  l!-lh,  a  similar  set  ot  amplitude  and  phase  delav  data  is  .InqiHr.ed  f-ir  a  p;rr  [)ior- 
quality  subarrays  at  NoltSAK  and  LAb.X,  It  and  i  t  n  speeim  Iv.  I  he.si  nlns  illii.s*!  'in  itii 
effect  of  incoherent  sc-nsors  on  the*  subai  i  av  iKanus.  hti*  lo  d.-st  i  ui  t r- >•  iim  rfi  n  n.  1  du  suh 
array  beams  are  (b’ficient  in  eer  fam  !•  eqiii  nc k-s,  i  ausnig  I’.e.aks  m  *'ii-  .-f'e.'ti’ai  i  aMo  .t  si  nsfo- 
to  beam.  Such  peaks  are  seen  betwemi  1  and  ?.  11/  in  this  figisre.  A  e  .rn  spoivimg  se  i?ti 0  in  t  ie 
phase  delay  data  is  also  seen. 

At  1  llz,  the  stand.anl  deviatifins  at  Hu-  I .  \  and  NfH.'sAl!  ;  uh.ai  1  ivs  art  latlu-r  similar  foi 

the  event  considered,  ‘lhe.se  dev i.a!  loris  are  stiown  m  I  ig  L-  17,  Ha  avi-r  .aL"  d.  vtafions  are 
2.4  dll  at  NOHSAlt  ami  1.9  dH  at  LAsA.  an  iii.sigm t lennt  diiii  ••  lue. 

If  we  examine  how  the  staiula rd  deviations  varv  vviHi  f,  eiiiu  ni  u  i  Ini  ifo-  .  whv  He  si  i;s- 
mograni  amplitudes  show  more  variation  at  \t)I{.SAK  than  at  I, ASA.  l  igure  il-ls  sli  nvs  Hu 
standard  deviations  for  toui-  NDIiSi  it  and  hmi  I.ASA  stiliarmvs.  Ilu  suL.ii  i’avs  fioni  hoih  at  car. s 
show  a  similar  increase  in  standatd  deviation  Irom  I  to  ?  H/.  an.J  no  clerii  diff.  lem  es  m  scan.-c 
can  t)c  observed  between  IjXSA  and  .NntiSAU.  IhuM  v«.r,  t  iv  nig  ai  NtdiSAU  vi-rv  ofii-u 

have  spectra  peaking  at  higher  frequeni  u-s  than  at  I  A.~,A.  Tlu  sc  signiLs  tlu  1  ,  hu  *  expern  rno- 
more  scatter  and  amplitudt*  variation  on  M)HSAH  records  than  .it  I  X.SA. 
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('.W.  I-'i  asier 


I-  U.Avl.  WII'l  iiM.I  I- \  I  I  l.irs^  .NOiJSAK 


!i.  "  •  |>i  I  I  \  I  l!ii‘  sff  iii.!  »  \  )!  short- pv-i  iod  F-wnvi*  amplitudes 

•ii  o)|,’,  \l>  V  IS  ill  tii  Mis?  ra  1 »  .  ‘  i  titf  uiifi  otaiit  <  uirftisioii  uas  iliat  iiace  auifd  ii  iirles  in  nnnos  ^ 
s<  »ini  •  ri  I  ■  '  '  ;  idt  inoasiimiatit  to  ade(|uateK  desi’iilie  tin*  amplitudt’  \  rtriat  ions  nefoss 

l.lt oris  all-  iiii'ti'i  .‘.a.  i  >  aliio  tii  tlif  ,  [M  i  tral  a iiipi 1 1 tide  va i  iat r'lis  at  NOIlSAIl  ns  a  function 
ol  .(/iiiiiitli  and  dist  itn  of  tin  iitconiitii;  i Muits.  I  In*  pui\.osr  ;d'  this  is  to  provide  a  data  base 
\sliit  li  i-aii  I  «  11,0  (1  III  1- on  lUMctioii  a  itli  itu  snl. array  time  ru'siduals  ohtninetl  In  Sheppard  (sei*  p.  12 
ill  lo  t  1  I  1  I  .It  It  rmitit  flit  rriisfal  striu  tuif  iiiidt  r  NOIvSAIt.  This  rtiuld  provide  an  indepentient 
.•hick  on  tilt  \()IJS\II  i.ootli  1  prop '.set!  it\  Katu  .-.t  rom  which  was  detliieed  from  studies  of 
lot  al  I  \  •  111  s  ill  ;iiid!;iav!a  . 

Mil  l  otiip'if  1 1  i  iiial  '  In  iii.  prc>i  iitcd  h*  M  i.s  to  forni  lK*nins  usitij^  the  center  sensors  of  each 
snhai  ra\.  aid  coeipii't  tin  apt  .  trnl  ratio  . .f  .  aeli  i  enter  sensor  to  tht'  oeam.  Kvents  were  chosen 
winch  all  ha',  t  rood  s  ir.rial  - 1  <  >- no  t  s,  ratios  at  a  i;oinnion  f:  t*(jiien(  y,  her«'  taken  to  be  1 

M'u.s  i..r.  I  li-.t  n  e'.t ms  lent-  In  t  ri  pr  'ccssed  tn  this  manner,  thev  art*  listed  in  Table  II- 1 
ariatir.'d  m  orili  t'  if  iiv  rt  asiiu'  hacK  ;r'inMith  troin  .\<>I{SAK  to  t*neh  t'piet'nter.  Figure  11-19 

show;  III.  t  oiiiir  111  at  loll  of  hat  k  a/innrths  for  these  eveius  on  a  map  fd  NOKSAlt. 

I  ir  I  rich  of  flu  ah oxt  t  ;i  n's.  In  si  <  ol  «  ai  h  subarrav  center*  sensor  were  tnpt*red  using  a 

tin- still  i!  eo.orn-  w  nr  I  -w  .  Me  tap'  lo  <!  signals  wire  then  l-ourit>r  transformed  and  sniot>thed 

rnrnn  in  d  ,  *  •  i-.o.,  jji  •.  n  nnam.  Mic  speetial  l  atro  o|  eneh  signal  to  its  corresponding  beam  was 
l•olllp!ttl  (!  at  1  It/  m  di  I  thel-,,  mn  phutt'd  on  a  map  of  NOHSAlt.  The  resulting  maps  were  autti- 
rnatii  all’,  t  otif oiireii  hv  .■oini  nti  r  "smg  a  srniple  scheme.  I  lie  arrnv  was  dividi'd  into  triangular 
tin  I  ts,  •  ill  h  iriauuh  Irr.  itig  thti  e  suharrav  posifnins  as  vertices.  'Mie  thret-  amplitudes  at  the 
*,  I  rfi  •*  .  V  '  rc  Hied  to  deftm  .i  plane  of  difming  amplitude  over  which  amplitiules  V’ary  linearly. 

I  ontotir  Inn  s  .irc  therelore  strairlit  lines  on  each  triangular  facet,  which  always  connect  with 
nher  .-oni  inr  lini  ni  atiia'  '•ill  frt  '  tn  I  In-  put  p-tse  »l  this  was  to  produce  simple,  unique  maps 
noi  snl'i*  I  !  '  I  iiiii  rpn  ’im  sm  lothing  tiv  tfic  analvsi. 

I  iL"  I  s  !  -  'ni  ii  inroiiL'h  it-  i  ari  the*  c  uit mred  m.aps  of  amplitudes  with  2-dH  contour  inter¬ 
val.-.  Il  l  ’  LOCI  -.n  III  't|i  o!.ji  •  ni-  ren-m.,  a/r.nitli  ot  eva n's  in  Talfh*  II- 1,  and  this  si— 
ipii  II*  •  IS  tin  in  c-t  iMi.tcsi'nc  'lie  I  •  l•^anllru  .  I*  lyures  1 1  -  2  hic  I  and  (d  )  a  re  a  mpl  itude  maps 

■  ir  iw  '  pi  I  sill!  I  I  I  \plo,;!  lit.  Ill  laasicrn  Kn/.fikh  and  ihev  are  quite  similar.  In  h’ig.  Il-20(h),  the 

ainplioidi  for-  ,  I’and.i  a  i  a ‘'tliipriki'  is  .shown,  lliis  ‘-veni  is  nearly  in  thi*  shadow  zone  at  a 
.ii.^i  iiiei  u  to*'  ,  VI  !  .1  pr  Mini  I  S  .in  .aniplitiidi-  innf)  with  man*  features  .similar  to  those  in 

i  .  ;.  li-.'in.  !  jino  'di  Me-  i  -  i  ans'  the  lianda  Si  a  I'.'ent  has  an  azimuth  elosi*  to  that  of  the 

Iv.i/.akt  I  VI  Id..  1  it'ufi  s  Il-2(iie.  .aitd  ifi  show  \  i  rv  similar  amplitude  patterns  for  two  Hindu  Kush 
*  ni  d  n'ia’.  I  S  V  ifh  fill!'  11  nt  di  pihs.  '  U*;irly  fti. -sc-  patti  rns  a t*e  caused  hv  crustal  and  mantle 
.fill,  iuii  nii'i-  I  NOl.VSAi.',  \it  the  pitti  riis  are  .so  sensitive  to  event  azimuth  that  Figs.lI-ZOfc) 

•iii-l  iM  :i  I  ■  .put.  .hpir’.'id  iioin  u  i  ati'i  nil,  the  difference  m  azimuth  In'tween  thi*  two  pairs  of 
.•V I  111  s  In  Ml  >  cihv  :,!i  Mil  2')  A  n.’nil.ii  lift  leulf  v  oi  •  ms  in  I  igs.  1 1- ZOfiM  and  f  i  1  which  show 
di  ..simil  .r  'i.itti  in."  s-n-  tw  -  •  \  .  nr  will  azimuths  differing  hv  onl*  10'. 

rti!..  nrfi  .'vts  iii.-d  s  iiie-  sm  .otIiniL'  dii  .a/tin.ith  inav  hi  required  in  order  t>  adequatclv 
lo  mIi  I  *.  ■  1  r  ii.-ne  !  on  I’l  .p  ins'  i ii  i hi  -i  fin  t  i.n  iiiuii  r  N<  )|{.SA It . 
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In  conjunction  with  this  work,  a  program  is  being  written  to  compute  the  three-dimensional 
diffracted  field  transmitted  by  realistic  models^ ^  of  the  MOIIO  under  NORSAR.  The  method  is 
based  on  an  application  of  Green's  theorem  described  by  Trorey^^  as  used  in  seismic  exploration 


work. 


C.  W.  Frasier 


F.  EFFECTS  OF  CRTSTAL  STIU’CITCRE  BENEATH  NORSAR 

ON  SHORT- PERIOD  TELESEISMS 

Finite  difference  solutions  for  teleseismiic  P-waves  emerging  into  an  irregular  boundary 
with  an  impedance  contrast  characteristic  of  the  Mohorovicic  discontinuity  beneath  NORSAR  have 
been  generated  to  determine  if  .such  a  structure  can  cause  the  amplitude  anomalies  observed  at 
that  array.^' ^ ^  The  shape  of  the  boundary  is  taken  from  Kanestrom.^^  His  model  is  based  on 
dT/dA  variations  across  the  array  and  agrees  with  Sheppard's  model^  based  on  the  travel¬ 
time  anomalies  across  the  array.  Depth  to  the  .Mohorovicic  discontinuity  for  the  model  is  con¬ 
toured  in  Fig.  11-21.  The  three-dimensional  structure  has  been  divided  into  2  two-dimensional 
structures  along  A- A'  and  B-B'  to  facilitate  computation.  Consequently,  the  wave  solutions 
will  not  be  valid  in  the  transition  regions. 

Along  A-A',  the  crust  thins  about  8  km,  in  isostatic  response  to  the  less-dense  rocks  com¬ 
prising  the  Oslo  Graben.^^  This  step- like  shape  is  shown  in  Fig.  11-22.  The  solid  line  is  taken 
directly  from  Fig.  11-21;  the  da.shed  line  is  the  approximation  used  in  the  computations.  Fig¬ 
ure  H-2  3  shows  the  velocities  and  densities  and  illustrates  the  initial  conditions  for  which  solu¬ 
tions  were  obtained.  Plane  harmonic  0.5- sec  P-waves  emerging  from  the  mantle  side,  with  angles 
of  incidence  ranging  from  normally  incident  to  iBO"  by  10*  increments,  were  propagated  through 
the  model.  Solutions  at  0.75,  i,  1.25,  and  1 .5  sec  for  normal  incidence  were  also  obtained.  The 
numerical  technique  used  to  solve  this  boundary  value  problem  is  explained  by  Landers.^^  The 
solution  is  approximate  in  that  backscattered  radiation  in  the  X-direction  is  not  computed.  The 
grid  sizes  and  the  numerical  mixing  parameter^^  D  are  also  given  in  Fig.  11-23.  Figure  11-24 
shows  an  equal-amplitude  plot  of  the  solution  for  normal  incidence  at  a  period  of  1.5  sec.  Each 
character  is  the  amplitude  of  the  wave  at  that  grid  point.  The  alphabet  has  teen  used  to  ex¬ 
tend  the  integer  s  over  a  range  of  1  to  36  (e.g.,  A  =  10).  Negative  amplitudes  are  plotted  as 
blanks  to  accentuate  wavefronts  (amplitude).  The  solution  at  38  km  above  the  step  (i.e.,  across 
NORSAR)  shows  the  offset  wavefront  (time-anomaly)  and  the  amplitude  variation  along  the  wave- 
front  caused  by  the  step-like  structure  below.  Figure  11-25  shows  the  time  and  amplitude  anom¬ 
aly  at  38  km  above  the  step  for  normally  incident  waves  with  different  periods.  l*he  time  delay 
(the  constant  phase  wavefront)  is  essentially  identical  for  each  period.  The  amplitude  anomaly 
at  1.5  sec  indicates  that  the  diffractions  from  the  sides  of  the  step  and  the  refraction  from  the 
step  are  producing  significant  amplitude  variations  at  38  km.  For  0.5-sec  waves,  the  amplitude 
is  that  of  a  smooth  shadow  zone  above  the  step  with  a  2-to-l  amplitude  anomaly.  Figure  11-26 
shows  the  amplitude  anomaly  for  0.5-sec  waves  emerging  at  various  angles.  With  minor  vari¬ 
ations,  the  overall  shape  of  the  curves  remains  the  same.  2-to-l  anomalies  with  sharp  transi¬ 
tion  zones.  Figure  11-27  shows  the  amplitude  of  the  S- waves  produced  by  conversion  at  the 
interface.  \^'hile  the  anomalies  are  very  distinctive,  these  waves  are  essentially  horizontally 
polarized  and  so  not  observable  across  the  array  of  vertical  seismometers. 

The  same  set  of  solutions  computed  for  the  North- South  model  wc.re  computed  for  the  Elast- 
West  model  shown  in  Fig.  11-28.  As  illustrated,  the  mantle  bows  up  into  the  crust,  producing 
an  anticlinal  structure.  Again,  the  model  is  consistent  with  isostatic  compensation  at  the  depth 
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of  the  MohorovMClc  discontinuity  for  tiic  Icss-densc  rocks  filling  the  Penman  (tslo  (Irahen.* 
Figure  11-29  gives  the  pertinent  numerical  value's  atid  illustrates  the  model  and  initial  eonditions 
for  which  solutions  were  ol.tained.  I•■ieur<•  11-  <0  shows  the  ef|iinl- ampl i tude  f)loi  for  a  nonnally 
incident  1 .5- sec-period  l*-\\ave.  The  shape  anfl  depth  of  tiu'  anticline  »k*focu.s  the  incident  wave 
producing  a  time  advanced  lower  amplitiuie  wave  over  ihe  i  i*nter  c)f  the  stnu.  t m  e,  Fiyui  e  II-  il 
shows  the  anomaly  as  function  of  period  and  the  time  eh  lav  which  is  essentiallv  the  same  at  all 
periods  Once  again,  at  longer  periods,  rapidlv  varving  amplitude  anomalit's  of  about  4  ?  mag¬ 
nitude  beco  e  a  W'cll-defined  pair  of  shadow  /ones  with  amjilitude  anomalies  of  about  2:  1  at 
short  periods.  I'igure  11-32  shows  the  variations  in  amplitudes  for  various  iticidenei'  angles. 

In  summary,  rather  simple  structures  at  the  mantle-erust  interface  can  c’aiisf'  2:1  ampli¬ 
tude  variations  in  teleseismic  P-waves  over  di.stances  as  small  is  10  km.  'Ihe  si/e  of  the 
anomalies  should  decrease  with  increasing  period;  however,  tl.v  variations  at  these  periods 
may  not  be  as  .smooth.  Hecause  of  the  rapid  variation  of  amplitude  over  small  distances,  am¬ 
plitude  anomalies  observed  at  the  relativelv  widelv  spared  suharravs  ai'ross  NOKSAl!  should  not 
be  expected  to  v’ary  smoothly  for  events  with  smoothly  <  hangmu  epii  eutral  eo- ordinal es.  A  de¬ 
tailed  model  of  the  crust  and  upper  mantle  htmeath  the  array  roidd  fit-ovide  deterministic  travel- 
time  corrections  but  W'ould  be  necessary  to  :ompute  waveform  e.orrertions  for  improving  the 
array  location  and  detection  capabilities. 

'r.  Ik  Landers 


G.  PHASE  VELOCITY  AND  AZlML'TIf  DEVIA'l  lONS  ACKOS.S  l.A.SA  AND  NOKSAlt 
FOa  LOCAL  EVENTS 


The  source  of  the  large  amplitude  and  phase  anomalies  acros.s  1..ASA  and  NOlLS.AIt  has  not 
yet  been  characterized  adequately.  Various  author  have  shown  that  the  effects  of  a  laterallv 
varying  sedimentary  section,^  ^  an  irregular  Mohorovicic  discontinuity/'^^  a  random  seattering 
mechanism  in  the  upper  mantle,  ’  and  a  laterally  varying  structure  in  the  L  wer  mantle*  can 
explain  to  a  certain  degree  the  observed  fluctuations  in  the  deta.  While  it  is  most  likelv  that  all 
of  the  above  contribute  to  the  overall  behavior  of  seismic  waves  ob.^erved  across  the  arrays,  it 
is  of  primary  interest  to  determine  which  ones  have  first-order  effeets  on  slownes.s  and  travel 
times  if  the  use  of  these  arrays  as  arravs  is  to  be  maximized,  lor  example,  the  imprf)vement 
in  detection  level  could  be  as  much  as  an  order  of  magnitude  if  station  eorreottons  were  (lefer- 
Tuinistic,  since  a  subarray  beam  at  NOltSAH  may  l)o  no  larger  than  the  seismogram  from  a 
single  element  of  that  subarray.  To  ai<l  in  determining  the  nature  of  the  crust  and  uppermo.st 
mantle  beneath  the  arrays,  a  study  of  local  and  near- regional  events  has  been  undertaken.  .Such 
data  have  a  distinct  advantage  over  teleseismic  data  in  that  wave  patlis  traverse  only  the  region 
under  study.  The  use  of  an  array  of  subarrays  also  has  the  advantage  f)f  having  local  values  of 
phase  velocity  and  azimuth  as  well  as  travel  time  in  the  data  base. 

Tables  11-2,  11-3.  and  11-4  show  data  for  a  local  and  a  near- regional  ev«  nt  at  I  ASA  and  for 
a  local  event  at  NOKSAK.  I  lie  phase  velocity  C  and  azimuth  A/.l  at  each  sutiarray  vv*tc  deter¬ 
mined  by  time  picking  the  outer  6  elements  at  NOH.S.AH)  and  center  instrument  and  least-.s.juares 
fitting  a  plane  wave  to  the  picks.  The  rms  time  error  At,  as  well  as  the  rms  error  in  phase 

velocity  AC  and  azimuth  AAZl  arc  al.so  given.  AC  and  AAZl  are  computed  from  the  relation- 
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ship  given  by  Kelly 


AAZl  = 


C 


AT^C 
N  VAH  X 
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TAIil.K  n-2 

at  IS.Q  N.  14:18:00 


1 

AAZl  (rms) 

I  AZI 

i  Al  (rnist 

AP  (rms) 

(deg) 

\ 

^  I  lUH.H 

'  (1.01 0't 

0.14 

1.3 

>  210.'^ 

!  0  0  403 

0.39 

3.4 

203.7 

i  0.00‘»4 

0.07 

0.6 

1H7.3 

!  0.0160 

0.11 

0.9 

P>0.2 

0.f»214 

0.1  3 

1.2 

)  Ptd.H 

0.0079 

0.05 

0.8 

^  20  3.7 

0.02  42 

0.1  8 

1.4 

)  pid.O 

0.0106 

0.12 

l.i 

i  1H2.2 

1  0.0110 

0.07 

0.6 

3  20H.h 

0.(3186 

0.G9 

0.8 

1  20K.4 

0.0166 

0.1 1 

1.0 

3  !  1H4.9 

0.0183 

0.14 

! 

'  1  1H3.2 

0.0121 

0.08 

1  0.7 

3  pi^.o 

0.0166 

0.19 

1  1 

1.3 

^  1  2  30.2 

0.0384 

! 

2.3 

)  1  203. h 

i 

O.0H80  ’ 

1  6.19  i 

1.9 

?  '  P>6.2 

*  0.0344 

!  0.2  3 

2.0 

3  204.9 

0.0201 

0.20 

1.8 

1  i  267.6 

j  0.0696 

i  0.70 

5.2 

1  '  110.2 

0.0270 

'  0.21 

1  1.8 

I  ^  172.3 

'  0.062  8 

!  0.40 

!  4.; 

whcrf*  .\  IS  Uh*  MumIxT  ')t  rlnnrnis  and  \  Al{  \  is  tfie  rnis  %'alue  of  one  of  the  orthogonal  com- 
poinaits  <if  tlu*  Stan  III  I  Ih**  derualion  assumos  a  normal  riistT  ihution  of  travel-time 

errors  and  a  svinmetni-  arrav.  f  or  the  LASA  eonftguratton  (hexagonr  1 

.X/.,  til.-' 
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and  for  N'DhSAH  (petilag  >na  1 », 
AA/d 


A  I  L 
(•  ’  l.4d  H 


ivhere  H  IS  the  dtstatn  e  fr  *iit  die  t  enter  to  an  elenteni  of  the  outer  ring.  'I'he  quantizing  error, 

+  see,  leails  to  a  A'l’  of  n.d|  nnd  a  nominal  AAZl  and  A(*  of  1  ’  an!  0.1  kni/scc,  respectively 

tor  both  ai  ra%s.  As  sec'ti  from  1  aides  11-2  through  *’  1,  the  arr(\als  often  differ  significantly 
f  1(011  plane  \va’.  es  hv  more  than  (piantizing  t:an  aeeount  for;  however,  in  gt  neral,  arrivals 
are  nearly  planar  whereas  P^  arrivals  art*  somewhat  w*»rse.  Ttgurv'  11-  3^  shows  the  azimuths 
of  the  K-  and  1- rings  at  I.ASA  for  the  event  m  I  n'n’e  11-2.  i  he  epicenter  is  well  defined  in  an 
area  rd  ahoiit  28  km.‘'  The  P.^,  arrival  at  I'l  is  ahotn  10'’  in  error.  Figure  11-34  shfnvs  the  azi¬ 
muths  for  th*  NOK.SAIl  event  given  in  Table  II- 4.  The  epicenter  here  is  not  nearly  as  well  de¬ 
fined,  however,  excluding  tlie  suharravs  Ivmg  to  the  south  and  east  of  the  dashed  line  in  the 


TABLE  11-3 

EVENT  AT  44.3“N.  105.1  “W.  11:44:20.4 


Subarray 

C 

AZI 

— 

AT  (rms) 

AC  (rms) 

AAZI  (rms) 
(deg) 

AO 

8.06 

171.3 

0.0166 

0.18 

1.2 

B1 

9.13 

169.3 

0.0119 

0.16 

1.1 

B2 

9.89 

171.6 

0.0435 

4.0 

B3 

7.22 

161.2 

0.0242 

0.21 

1.7 

B4 

7.27 

169.3 

0.0395 

0.34 

2.8 

Cl 

7.54 

172.6 

0.0269 

0.25 

1.9 

C2 

7.75 

170.2 

0.0386 

0.38 

2.9 

C3 

9.75 

185.4 

0.0389 

0.61 

3.5 

C4 

7.96 

173.0 

0.0097 

0.7 

D1 

7.10 

179.1 

0.0240 

0.20 

1.6 

D2 

7.89 

163.7 

0.0233 

0.24 

1.8 

D3 

8.32 

172.3 

0.0513 

0.59 

4.0 

D4 

6.86 

168.8 

0.25 

2.0 

El 

8.22 

173.9 

0.0207 

0.23 

1.5 

E2 

8.62 

163.9 

0.0065 

0.08 

0.6 

E3 

9.27 

161.1 

0.0274 

o.i*- 

0.8 

E4 

11.40 

164.0 

0.0275 

O.f  9 

2.9 

FI 

7.93 

175.1 

0.0339 

0.35 

2.5 

F2 

7.81 

166.7 

0.06 

0.5 

F3 

9.22 

1 46.6 

0.0354 

3.0 

F4 

7.57 

166.8 

0.0256 

0.24 

1.7 

lABLE  n-4 

EVENT  AT  59.9 ‘N.  H.4*E,  16:44 


Subari'ay 

C 

AZI 

AT  (rms) 

AC  (rms) 

AAZI  (rms) 
(deg) 

01 A 

6.69 

165.9 

0.0036 

0.25 

0.22 

OIB 

6.70 

166.6 

0.0171 

0.12 

1.04 

02  B 

5.96 

174.5 

0.0297 

0.16 

i.60 

03B 

6.11 

175.3 

0.0102 

0.06 

0.56 

04B 

6.23 

173.8 

0.0187 

0.11 

1.05 

05B 

6.73 

160.4 

0.0158 

0.11 

0.96 

06B 

6.74 

149.9 

0.0206 

0.15 

1.26 

07B 

6.72 

160.2 

0.0334 

0.24 

2.03 

02C 

6.81 

169.3 

0.0240 

0.16 

1.48 

03C 

6.60 

185.1 

0.0087 

0.06 

0.52 

mm 

6.47 

188.8 

0.0087 

0.06 

0.51 

BH 

6.11 

189.4 

0.0176 

0.10 

0.97 

mm 

6.34 

185.3 

0.0152 

0.10 

0.87 

07C 

6.58 

184.2 

0.0192 

0.13 

1.14 

08C 

6.07 

157.6 

0.0331 

0.19 

1.82 

09C 

6.82 

142.3 

0.0300 

0.22 

1.85 

IOC 

6.73 

137.1 

0.0320 

0.23 

1.95 

lie 

8.14 

142.0 

0.0350 

0.37 

2.58 

12C 

7.30 

139.9 

0.0329 

0.28 

2.17 

13C 

6.75 

148.2 

0.0300 

0.22 

1.83 

figure  (2B,  3B,  5C,  4B,  5C,  8C,  and  7C),  the  epicenter  can  be  located  in  an  area  of  about  iOO  km.^ 
All  the  aforementioned  subarrays  then  have  large  easterly  azimuth  anomalies.  Figure  11-35 
shows  the  phase  velocity  measurements  for  the  same  event.  The  circled  values  are  inose  for 
the  subarrays  with  the  azimuth  anomalies.  They  are  anomalously  low  in  phase  velocity.  The 
upper  crust  beneath  NORSAR  seems  to  be  divided  into  at  least  two  provinces.  If  the  difference 
in  crustal  velocities  (6.1  km/sec  average  for  the  SE,  6.7  km/sec  otherwise)  persisted  over  30  km 
in  depth,  a  travel-time  anomaly  of  0.4  sec  would  be  expected.  In  Fig.  11-36,  the  phase  velocity 
at  the  different  subarrays  for  the  two  LASA  events  is  plotted  with  error  bars  of  ±AC.  The  var¬ 
iation  in  Pjj  values  is  considerably  larger  than  for  values.  Whereas  at  NORSAR  the  azimuth 
and  phase- velocity  variations  indicate  significant  lateral  variation  in  the  crust,  the  values  at 
LASA  indicate  that  the  greater  variation  in  structure  lies  wear  the  Mohorovicic  discontinuity. 

Continuing  analysis  of  the  large  number  of  local  strip- mining  explosion  records  at  LASA 
and  near  events  at  NORSAR  should  lead  to  a  better  understanding  of  the  ^rust  and  mantle  be¬ 
neath  these  arrays,  and  that,  in  turn,  may  lead  to  improved  array  capability. 

T.  E.  Landers 
H.  Bungum 
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H.  LOCATING  ABILITY  OF  LASA  AND  NORSAR 

Large  seismic  arrays  such  as  Luf'  SA  and  NORSAR  have  been  used  singly  to  locate  seismic 
events  on  a  routine  basis  for  several  y^ars,  but  few  attempts  have  been  made  to  use  two  or 
more  arrays  simultaneously  for  this  purpose.  The  event  location  algorithm  described  in  the 
last  SATS^  enables  one  to  calculate  the  theoretical  resolving  power  of  any  network  of  arrays 
and/or  single  stations,  and  we  present  here  as  an  illustration  of  this  ability  the  confidence  re¬ 
gions  expected  when  LAS^  and  NORSAR  are  used  separately  and  in  conjunction  to  locate  events 

19 

in  central  Asia.  Gjoystdal,  et  analyzed  this  problem  using  a  Monte  Carlo  technique  to 
simulate  random  errors. 

Figure  11-37  shows  the  theoretical  resolving  power  of  LASA  for  shallow  events  in  central 
Asia.  These  calculations  are  based  on  the  assumption  that  LASA  can  measure  slowness  with 
a  standard  deviation  of  0.05  sec/deg,  and  thus  represent  the  "best  possible”  situation  in  which 
virtually  perfect  site  corrections  for  *hz  array  sensors  are  available.  Without  such  corrections, 
systematic  location  errox  s  would  of  course  be  much  larger  than  those  illustrated.  The  regions 
shown  are  formal  50-percent  confidence  ellipses  for  a  normal  distribution  of  observational 
errors.  Regions  of  greater  confidence  have  deliberately  been  avoided,  since  their  size  would 
depend  critically  on  the  details  of  the  "tail”  of  the  distribution  of  measurement  errors.  The 
numbers  given  to  the  right  of  some  of  the  ellipses  give,  respectively,  their  aspect  ratios  and 
areas  (in  thousands  of  km^).  Several  points  about  Fig.  11-37  should  be  noted.  The  loss  of  res¬ 
olution  near  the  bottom  of  the  figure  as  the  core  shadow  is  approached  is  obvious.  The  over¬ 
lap  of  regions  of  radically  different  size  indicates  that  the  linearization  of  the  equations  implicit 
in  these  calculations  is  not  a  valid  approximation  in  the  large  ellipses.  The  large  differences 
in  size  between  some  adjacent  ellipses  are  caused  by  irregularities  in  the  theoretical  slowness- 
vs- distance  curve,  and  probably  are  not  meaningful. 

Figure  11-38  shows  the  corresponding  confidence  regions  for  NORSAR  (assumed  to  have  a 
resolution  of  0.1  sec/deg).  Some  of  the  apparent  difference  between  the  ellipticities  of  the  re¬ 
gions  on  Figs.  11-37  and  11-38  are  the  result  of  distortion  introduced  by  the  map  projection  (cf. 
the  regions  at  65*N,  140’E).  The  regions  in  the  upper-left  corner  of  Fig.  11-38  should  not  be 
taken  seriously,  since  they  are  less  than  25*  from  NORSAR,  w'here  the  travel- time  curve  is 
multi-branched  and  strongly  varying  regionally. 

Figures  11-39  and  11-40  show  the  confidence  regions  expected  if  LASA  and  NORSAR  are  used 
together  to  locate  events.  Figure  11-39  is  based  on  an  uncertainty  in  absolute  travel  times  of 

I. 5  sec,  which  probably  corresponds  fairly  well  to  the  accuracy  obtainable  with  simple  tables. 

In  Fig.  11-40,  this  figure  has  been  reduced  to  0.5  sec,  which  might  be  obtainable  if  regional 
variations  were  included  in  the  tables.  It  is  evident  that  accurate  arrival- time  data  are  poten¬ 
tially  much  more  useful  than  slowness  measurements,  although  the  earth  must  be  "calibrated” 
to  fully  take  advantage  of  them.  It  is  also  clear  from  Fig.  11-40  that  the  optimum  location  for 
another  low-threshold  station  would  be  somewhere  in  southern  Asia. 

B.  R.  Julian 
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crosscorrelation  coefficient 

Fig. 11-11.  Theoretical  relationship  between  amplitude- relative- to- phase 
fluctuation  power  ratio  and  crosscorrelation  coefficient  between  amplitude 
and  phase  fluctuation,  and  the  90- percent  confidence  limits  of  measurement. 
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Fig.  11-12.  Comparison  of  estimate  with  theoretical  curve  for  transverse 
autocori  elation  coefficient  of  amplitude  fluctuations. 


Fig.  11-13.  Compari.son  of  estimate  with  theoretical  curve  for  transverse 
autocorrelation  coefficient  of  phase  fluctuations. 
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Fig.  11-14.  Typical  subarrays  at  NORSAR  and  LASA.  All  sensors  of  2C 
at  NORSAR  are  shown,  whereas  for  F3  at  LASA  only  outer  ring  and  cen¬ 
ter  sensor  are  shown. 
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Fig.  11-15.  Amplitude  and  phase  delay  statistics  of  spectral  ratio  of  single  sensor 
to  subarray  beam  using  presumed  explosion  in  Eastern  Kazakh,  6  June  1971.  Solid 
lines  indicate  mean  amplitude  and  phase  delay,  and  dashed  lines  indicate  .scatter 
within  a  standard  deviation  of  mean.  Subarrays  6('  at  NOHSAR  and  F3  at  LASA 
are  shown. 


NOHSAR  1C  subarray  LASA  C*  SUBARRAr 

Fig.  II-  16.  Aniplltude  and  phase  delay  statistics  as  described  for  Fig.  II-  15. 
In  this  case,  subarrays  1C  at  NOHSAR  and  C4  at  LASA  arc  shown. 
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NORSAR 


Fig.  II-20(a-k).  Spx‘»ctral  ratio  at  1  Hz  of  subarray  center  sensor  to  array  beam 
at  NORSAR  contour  interval  is  2  dB,  using  events  listed  in  Table  II- 1. 
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Fig.  Il-20{a-k).  Continued. 
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Fig.  11-23.  Velocities  and  initial  conditions  for  which 
step- model  solutions  were  generated. 
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Fig.  11-24.  Kqual-amplitude  plot  of  solution  for  1.5-sec  period  plane  P-wave 
normally  incident  on  North-South  mt>del. 
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Fig.  II-2S.  Amplitude  and  time-delay  Fig.  11-26.  Amplitude  anomah.'»s  for  step 
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Fig.  11-27.  S-wavc  amplitudes  for  North-South  model  for  P-waves  incident 
at  0"  and  *30'. 


Fig.  11-28.  Cross  section  along  B-B' 
and  approximating  anticline  model. 
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Fig.  11-29.  Velocities  and  initial  conditions  for  which  anticline 
model  solutions  were  generated. 
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Fig.  11-30.  Ek)ual- amplitude  plot  of  solution  for  1.5-sec  period  plane 
P-v/ave  normally  incident  on  East- West  model. 
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P-WAVC  AMPtITUOES 


Kig.  11-  3^.  Azimuths  of  first  arrivals 
at  IC-  and  K-ring  LASA  subarrays  for 
event  in  Table  II-2. 
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Fig.  11-  34.  Azimuths  of  first  rrrivals  at  C-ring  NORSAR 
subarrays  for  event  in  Table  11-4. 
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Fig.  11-37.  Theoretical  50-percent  confidence  ellipses  for  LASA, 
based  on  assumed  resolution  of  0.05  sec /deg.  Numbers  are,  re¬ 
spectively,  aspect  ratios  and  areas  (thousands  of  km^)  of  ellipses. 
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Fig.  11-38.  Theoretical  50- percent  confidence  ellipses  of  NORSAR, 
based  on  an  assumed  resolution  of  0.1  sec/deg. 
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Fig.  11-39.  Same  as  Figs.  11-37  and  11-38,  for  LASA  and  NORSAR 
used  together,  with  travel- time  curve  accurate  to  1.5  sec. 
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Fig.  11-40.  Same  as  Fig.  II- 39,  with  travel-time  curve 
accurate  to  0.5  sec. 
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UI.  SKISMK'  DISCIUMINATION 


A  •  «:ti:(  tk)n  ok  intkkkfmiing  itAvi.Kir.ii  waves  at  noksak 

iv  is  well  known  that  a  powerful  technique  for  discriminating  between  natural  events  and 
undeigi‘r)un(l  niu  lcar  explosions  is  the  M^-nij^  method,  based  on  the  relationship  between  the 
surface-wave  ^^nd  the  body-wave  magnitude.  Considerable  effort  has  been  expended 
on  the  problem  of  magnitude  definitions,  and  especially  on  the  definition  of  yielding  the  best 
possible  separation  between  earthquakes  and  underground  explosions  on  an  M^-m^  diagram, 
cf.  Kef.  1.  However,  with  the  present  dcpl(\yment  of  sei.'^mic  observatories  and  arrays,  the 
main  limitation  of  the  method  is  (piite  often  that  no  Kayleigh  waves  can  be  found  for  an  event 
which  is  detected  and  located  on  the  basis  of  short-period  P-wave  recordings.  One  of  the  rea¬ 
sons  for  this  lack  of  a  detection  of  surface  waves  is  that  their  energy  level  may  be  below  that 
of  the  ambient  background  microseismic  noise.  In  this  case  the  discriminant  may  still 

be  employed,  above  a  certain  body-wave  magnitude  level,  provided  that  discrimination  can  be 
based  on  negative  evidence. 

However,  quite  often  there  is  a  serious  detection  problem  caused  oy  the  long  persistence 
of  Kayleigh-wave  coiia  of  other  events,  as  pointe'^’  out  by  ('apon  and  Evernden.^  This  situation 
is  complicated  by  the  fact  that  rapid  lateral  variations  in  phase  velocity  frequently  cause  severe 
multipath  propagation  problems,  as  pointed  out  by  Capon.^  The  purpose  of  the  present  study 
was  to  investigate  the  amount  and  the  characteristics  of  such  multipath  propagation  for  Rayleigh 
waves  recorded  at  NOK.SAH,  and  to  Cf)nsi<ler  the  problem  of  detecting  one  Rayleigh  wave  in  the 
presence  of  the  coda  of  another  Rayleigh  wave.  The  long-period  (I.P)  array  at  NORSAR  consists 
of  a  three-component  set  of  seismometers,  located  within  eai'h  of  the  22  subarrays,  with  peak 
respon.se  at  about  25  sec.  The  aperture  of  NOK.SAR  is  approximate!}'  110  km. 

The  main  analysis  technique  used  in  this  work  was  the  high- resolution  (HR)  frequency- 

A 

wavenumber  spectrum  analysis  method  describeil  in  detail  by  Capon.  The  technique,  which 
essentially  reveals  information  about  the  distribution  of  power  as  a  function  of  frequency  and 
wavenumber,  has  been  applied  to  successive  nonoverlapping  200-sec  intervals,  at  periods  of 
40  and  20  sec,  starting  at  the  onset  time  of  the  lO-sec -period  Ktiyleigh-wave  group.  The  anal¬ 
ysis  of  such  short  time  intervals  would  normally  entail  a  serious  windowing  effect,  causing 
energy  fiom  frequency  bands  other  than  those  under  analysis  to  affect  the  results,  cf.  Capon. ^ 

In  order  to  avoid  this,  a  prefiltoring  operation  has  been  applied  to  the  data  before  the  HR  anal¬ 
ysis.  The  filter  used  in  the  40~sec  period  analysis  was  a  fifth-order  Chebyshev  low-pass  filter 
with  3-dB  response  at  35  sec,  while  for  the  20-sec  analysis  a  third-order  Butte^worth  bandpass 
filter  was  employed  with  3-flB  responses  at  25  and  16.7  sec.  Both  filters  are  rolling  off  very 
sharply,  thus  virtually  eliminating  all  problems  connected  with  frequency  windowing. 

The  HR  method  was  applied  to  prefiltered  data  from  the  I.P  vertical  array  at  NOR.SAR  for 
15  events,  each  with  >  5.2,  located  at  various  azimuths  and  distances  from  NOR.S.AR.  There 
were  no  events  occurring  close  in  time  to  the  cho.sen  events,  according  to  the  l  SCGS  monthly 
summar  reports  and  the  NOR.S.AR  bulletin. 

The  method  proposed  by  Capon  and  Eivernden^  for  the  detection  of  interfering  Rayleigh 
waves  at  LA.SA  has  been  considered  also  for  NOR.SAR.  This  method  entails  the  application  of 
the  HR  method  to  the  prefiltered  data  from  the  I.P  vertical  array  at  LA.SA,  to  determine  the 
angle  of  arrival  of  the  fundamental-mode  Itayleigh  wave  at  the  appropriate  arrival  time.  If 
this  angle  and  the  <  rrival  time  are  reasonably  close  to  the  expected  values,  then  a  detection  is 


said  to  occur.  In  order  to  evaluate  the  effectiveness  of  this  detection  method,  it  is  necessary 
to  measure  the  rate  of  decay  of  the  Rayleigh-wave  coda  power  level  at  NORSAR.  These  results 
are  given  in  Figs.  Ill-l(a)  and  III-2(a)  for  the  20-  and  40-sec-period  groups,  respectively.  The 
data  in  these  figures  were  used  to  determine  the  minimum,  average,  and  maximum  detection 
thresholds  for  the  detection  method  as  a  function  of  the  arrival  time  of  the  surface  waves  of  the 
sought  event,  and  the  results  are  given  in  Figs.  Ill-l(b)  and  III-2(b)  for  the  20-  and  40-sec -period 
groups,  respectively.  In  obtaining  these  figures,  detection  thresholds  of -12  and -6  dB  below 
the  coda  power  level  were  used,  respectively,  as  determined  experimentally  for  the  20-  and 
40-sec-period  group  Rayleigh  waves  at  NORSAR. 

It  is  seen  from  Figs.  Ill-l(a)  and  III-2(a)  that  after  about  7  minutes,  the  40-sec-period  group 
power  has  dropped  about  8  dB  more  than  that  for  the  20-sec  period.  These  results  should  nor¬ 
mally  mean  that  by  analyzing  at  a  40-sec  rather  than  at  a  20-sec  period,  the  detection  of  the 
Rayleigh  waves  would  not  be  obscured  as  often  by  an  interfering  event,  provided  there  is  equal 
signal  energy  at  these  two  periods.  Another  argument  that  also  favors  the  40-sec-period  group 
is  the  observation  that  there  is  more  severe  multipath  propagation  at  '''»n  the  other  side, 

however,  is  the  fact  that  the  HR  method  performs  better  at  20  sec,  due  to  the  ui.,  )ved  array 
response  at  this  period.  In  addition  the  signal-to-intcrference  ratio  may  actually  be  better  at 
20  sec,  if  there  is  enough  signal  energy  to  compensate  for  the  increased  amount  of  interfering 
event  energy  at  this  period,  relative  to  40  sec.  A  general  conclusion  as  to  which  frequency  to 
use  therefore  cannot  be  reached,  and  the  detection  at  these  two  periods  should  be  considered 
simultaneously. 

The  HR  analysis  of  the  15  events  provides  information  concerning  the  azimuthal  deviation 
from  the  great  circle  path  for  all  time  frames  in  which  the  phase  velocity  measured  can  be 
reasonably  associated  with  fundamental-mode  Rayleigh  waves.  This  acceptable  range  of  phase 
velocities  was  3. 4-4.0  and  3. 5-4. 2  km/sec  for  20-  and  40-sec-period  Rayleigh  waves,  respec¬ 
tively.  These  data  concerning  angles  of  arrival  and  group  delays  can  be  used  to  infer  the  actual 
propagation  paths  for  the  multipath  propagation  of  Rayleigh  waves,  cf.  Capon.^  The  results  for 
the  propagation  paths,  at  both  20-  and  40-sec  periods,  for  the  7  December  1971  Bouvet  Islands 
and  18  December  1972  Zambia  events  are  shown  in  Figs.  HI- 3(a)  and  (b)  and  III-4(a)  and  (b', 
respectively.  These  data  seem  to  agree  with  the  corresponding  results  for  LASA  given  by 
Capon,^  since  it  appears  that  in  many  cases  the  multipath  propagation  can  be  associated  with 
refractions  and  reflections  at  continental  margins,  as  well  as  scattering  from  small-scale 
tectonic  regions  of  the  earth. 

It  is  now  shown  how  the  proposed  detection  method  has  been  used  at  NORSAR  to  detect  inter¬ 
fering  Rayleigh  waves.  A  number  of  independent  sources  reported  that  an  atmospheric  nuclear 
explosion  occurred  at  the  Lop  Nor  test  site  in  China  on  27  .Tune  1973.  This  explosion  has  been 
identified  by  the  USCGS,  on  the  basis  of  seismic  observations,  as  the  27  June  1973  Southern 
Sinkiang  event,  and  has  been  used  as  a  good  example  of  the  performance  of  the  proposed  detec¬ 
tion  method.  The  reason  for  this  is  that  the  event  occurred  during  the  time  that  an  earthquake 
swarm  was  taking  place  in  the  Kurile  Islands,  Hokkaido  region,  and  the  expected  arrival  time 
for  the  explosion  Rayleigh  waves  coincided  exactly  with  the  Rayleigh  waves  from  one  of  the 
earthquakes  in  the  swarm,  namely  the  27  June  1973  Hokkaido,  Japan  event.  In  addition,  an 
event  with  epicentral  parameters  which  were  quite  close  to  those  for  this  interfering  event,  and 
which  occurred  about  30  minutes  before  it,  was  chosen  ^s  an  ideal  reference,  namely  the 
27  June  1973  Kurile  Islands  event. 
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In  ucrordanrn  with  the  pi-ocodii i*c  outlincti  previously,  both  the  40-  and  ZO-sec-period  groups 
have  boon  analyzed  using  the  llll  method,  as  applied  to  the  pretilLered  data.  The  results  for  the 
40-sec'-peri()d  group  are  shown  in  Fig.  IIl-S(a-<n  where  thf?  UK  specti-a  from  four  consecutive 
liOO-sec  lime  intei  vals  are  deouded.  The  IIK  analvi^is  was  started  about  3  minutes  before  the 
computed  ari  i\  al  time  for  the  40-see-perio(l  Kay leigh-wave  group  for  the  desired  event,  i.e., 
the  27  .lime  197  3  .Southern  Sinkiang  c'vent  (Chinese  Atmospheric  .\uclear  Explosion).  The  main 
40-.sec-period  group  arrives  duiang  tne  secf)nd  time  frame,  depicted  in  Fig.  ni-5(b),  and  is  also 
present  in  the  third  time  frame,  given  in  f'ig.  !!!-S(c).  befoie  it  fades  aw'ay.  There  is  no  multi- 
path  propagation  obset'ved  in  these  two  time  frames,  since  the  angle  of  arrival  coincides  almost 
precisely  with  the  azimuthal  diiection  toward  the  flesired  event.  It  can  be  concluded,  thereiore, 
that  the  40-.sec-perio{l  Kayleigh-wave  group  of  the  flesired  event  has  been  detected,  since  it  has 
been  shown  to  c.*"*'ive  at  NOK.SAK  at  the  appropriate  time  and  with  the  expected  angle  of  arrival. 

The  situ;’Jf)n  ff)r  the  20-sec-perifHl  gj*oup  is  much  more  complicated.  The  results  of  the 
UK  analysis  ff)r  the  reference  event,  aiul  for  the  desired  plus  interfering  events,  are  shown  in 
Figs.  Ill-6(a-fl)  and  in“7(a-d),  respectively,  and  the  twi>  figures  should  be  studied  simultaneously. 
The  UK  analyses  given  in  Figs.  III-6  and  flI-7  are  started  at  the  e.xpected  onset  times  of  the 
20-.see-period  group  of  the  reference  an<l  interfei  ing  events,  respectively.  In  Figs.  lll-6(a)  and 
III-7(a),  the  angle  of  arrival  agrees  well  with  the  azimuth  f)f  bf)th  the  reference  and  interfering 
events,  namely  abfuit  33".  The  secoufl  frame,  given  in  Figs.  lll-6(b)  and  lll-7(b),  shows  that 
the  angle  of  arrival  decreases  to  about  IH’  for  both  events,  possibly  inflicating  a  reHection 
somew'.*,ei'e  in  northern  .SiluM  ia,  and  the  surd'ace  waves  for  the  explosion  ar’O  lurt  seen  as  yet. 

In  Fig.  Ill-6(c)  there  is  an  arlditional  multipath  arrival  (at  about  dS’l  foi’  the  reference  event, 
while  in  Fig.  l!I-7(c'  fhe  main  energy  now  arrives  at  an  azimuth  *  espernding  to  tlie  direction 
towarfl  Lop  \f»r.  Since  this  time  agr-ees  will  with  the  expf'ctcfl  .  '’rival  time  for  the  20-sec- 
period  Kayleigh-wave  group  of  the  desircfl  event,  the  detection  of  the  surface  waves  for  this 
event  has  been  firmly  established.  The  results  in  Fig.  in-7(d)  indicate  some  energy  which  is 
probably  due  to  the  desired  event,  while  the  flata  in  Fig.  lll-fjfd)  show  only  some  weak  multipath 
component. 

In  this  case,  the  detection  of  the  Kayleigh  wave  of  the  desired  event  in  the  coda  of  an  inter¬ 
fering  event  has  been  accomplished  at  both  20  and  40  sec.  However,  it  is  clear  that  the  detec¬ 
tion  at  a  40-sec  period,  at  least  in  this  exanipie,  was  much  clearer  than  that  at  a  20-sec  period. 

J.  Capon 

II.  Hungum 

B.  PRECUKSOK.S  TO  pF  FKOM  DEEP  SOI  TH  AMERICAN  FAKTIIQI  AKES 

Occasionally,  when  processing  the  early  short-period  arrivals  from  deep  South  American 
earthquakes,  automatic  detection  schemes  used  by  the  large  arrays  do  not  correctly  identify 
the  pP  depth  phase.  These  schemes  interpret  an  intermediate  phase  between  P  and  pP  as  pP, 
and  identify  the  true  pP  as  a  second  event  from  the  same  region.  If  these  intermediate  phases 
are  of  a  deterministic  origin,  there  are  two  plausible  explanations  for  that  origin.  Eithe»‘  they 
are  the  P-waves  from  small  earthquakes  which  follow  the  larger  deep  events,  or  they  are  reflec¬ 
tions  off  some  boundary  between  the  focus  of  a  single  earthquake  and  the  surface  of  the  earth 
above  the  focus.  We  tend  to  reject  the  suggestion  that  these  phases  are  due  to  reflections  trom 
contrasts  beneath  the  recording  station  since,  in  that  case,  these  phases  would  be  generally 
observed  independent  of  the  geographic  region  or  depth  of  the  source.  We  shall  attempt  to 
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intcrpict  these  mtiTnuMitate  phases  hotwet‘!i  I*  an<l  pi*  as  reflections  i  'om  a  physical  contrast 
above  the  tle{’p  sonia  {'s  and,  following'  \\!iit<'oinb  and  ViKlerson,*’  introduce  the  nomenclature 
pdl’  where  d  is  the  depth  Ixuieatli  the  surface  of  that  contrast. 

I'imtre  lll-K(a)  shows  th.c  slio rt -period  sul)ai-ray  sums  «»f  the  three  outer  rings  of  l.ASA  and 
the  initial  arri\als  from  »n  (‘arflupiakt  locateil  S  J7  kir.  h  .leath  northern  Argentina.  .-\s  indicated 
on  the  figuri',  there  is  a  relatively  small  yet  < dear  an<l  <  «»ht  ent  arriv!  across  the  entire  array 
about  midway  between  i'  and  pH  Kigu.e  lII-H(h)  is  sim:  •  to  III-  .(at,  si  ).ving  array  traces 
for  an  earthijuake  about  KOO  km  south  of  that  in  Hl-H(al  at  a  depth  of  about  SHO  km.  Here,  aHhough 
the  phase  rnidwav  between  H  and  pH  exists,  it  is  not  as  clear  as  the  previous  example.  How'ever, 
there  is  a  larger,  intermediate  seconilary  arrival  ab  )Ut  30  sec  before  pH  in  thi.s  ca.ae  In  order 
to  verify  the  existence  of  thest*  pha.ses,  we  have  searched  t‘te  short -period  vertical  seismograms 
of  the  World  W  ide  .Standard  Seismog»*apliic  Network  (WW.S.SNi  of  the  two  events  in  Figs. 
and  (b)  and  four  additional  deep  carthipiakes  HI  ^  500  km)  under  nfirthwestern  Argentina.  .Mthough 
these  intermediate  phastws  cannot  be  rea<l  from  single  station  seismograms  as  confidently  as  they 
c.an  trom  array  data,  we  have  attempted  tfi  time  an  arrival  within  30  sec  before  pH  and  thela.;,est 
arrival  between  H  and  pH.  yet  not  part  of  the  H  coda,  ('are  was  taken  to  exclude  HcH  from  the 
study. 

'I'he  lesults  of  these  re.idings  .ire  shown  in  Fig.  111-9  where  the  time  differences  pP-pdP  of 
the  largest  ph.ise  (open  cirtdes)  and  phases  within  30  sec  c)f  pP  (black  dots)  are  plotted  vs  epi- 
ccrtral  di  »  inro,  'Fhe  s«di«l  lines  show  how  these  time  differences  should  vary  with  distance 
given  h  i/ontal  rellectors  at  various  depths  <l  above  the  source.  I'roni  Fig  111-9,  the  strongest 
pflH  p  -e  would  appear  to  come  from  a  boundaiw  between  200  and  300  krn,  say  250  km.  The 
Solid  Cots  sc’atter  abon?  the  tOO-km  depth  line  (they  were  picked  to  do  so)  and  do  not  show'  evidence 
fd’  any  well-defined  horizontal  reflector  at  that  ilepth.  'I'he  location  cf  the  eaiihcjuakes  used  ex¬ 
tends  south  a  long  longitude  f,  3  W  ft  «>m  22 to  2H‘.S  .and  isasscu  iated  with  a  rather  complicated 
tectonic  system  which  is  proh.ibly  not  wcll-rcpresentc«l  ity  plane,  horizontal  boundaries.  These 
qualification.,  may  account  for  some  of  the  scatter  observed  in  Fip.  Ili-^' 

There  is  one  ii  tercstinc  cnrrclati*)n  that  c.an  be  drawn  "'ctween  the  pha.se  p250H  and  the  re¬ 
gional  seismicity  pattern.  In  Fig.  1)1-10,  we  have  |>rojectcd  all  *.'r  r'hqual.es  of  rn^^  >  5.0  that 
oecurred  from  19f>7  to  t97v)  between  iO'S  and  30’.S  and  f»2’\\  and  7a'\\  onto  a  plane  of  constant 
latitude.  This  craiss  section  show  s  the  well-known  west-east  dip  of  Andean  seismicity  and  the 
gap  in  this  activity  between  300-  and  500-kni  <lepth.  It  is  well -establisbed  that  this  general 
.seismicity  patlei  n  i.s  as..f<>ciated  with  the  plunging  of  the  riceanic  N’az.ca  plate  beneath  the  .South 
AmcriCciii  continent.  Isacks  and  Molnar^’  liavc  discussed  seismicity  gaps  as  seen  in  Fig.  111-10 
in  terms  of  physical  gaps  in  the  desceiuiing  litho.spherc,  and  as  representing  stress  changes  in 
the  flab  fr' m  extension  to  ronipres..,ion  as  a  ftincti<*ii  of  uepth.  The  riata  of  Figs.  lll-H  and  111-9 
indicate  that  a  rellei  ting  iKuind.arv  may  exist  under  we.stern  .South  .Vi  icri-^a  at  a  depth  that  coin¬ 
cides  with  the  beginning  of  the  larg*  gap  in  a  rathe,  narrow  seismicity  zone.  .Since  this  boundary 
must  be  rather  w  nlcspread  in  ordei*  tf)  reflect  waves  from  earth(|uakes  hOO  krn  apar.,  its  exist¬ 
ence  might  'ndicate  that  i-ausi?  of  the  stnkinglv  dis^'oiitininais  sei  .mirity  zone  is  not  of  the  de¬ 
scending  plate  it'.clf  hut  of  the  mantle  into  which  it  intrudes. 


.1.  IF  Fil.son 
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C.  COMPARISON  OF  pP-P  DELAY  TIMES  FROM  SI' RFAC E-ZERO 
AND  TELESEISMIC  RECORDS  OF  NUCLEAR  EXPLOSIONS 


Many  attempts^"^^  have  been  made  to  estimate  th'.  pP-P  delay  times  for  nuclear  oxplosif)ns 

using  teleseismic  short-period  seismograms,  and  from  these  delay  times  to  estimate  depth  of 

burial.  One  technique  is  to  look  for  repeating  nulls  in  the  spectrum  of  the  P-wa/e  train  and 

interpret  them  as  a  destructive  interference  between  P  and  pP  phases.^*®  Another  method  is  to 

deconvolve  the  combined  effect  of  seismometer  and  eartl  attenuation  from  the  seismogram  in 

the  time  domain,  and  obtain  the  pP-P  delays  by  inspection  of  the  deconvolved  rec  rds.*^’^® 

Recently,  homomorphic  deconvolution  has  been  utilized  to  extract  pP,  P,  and  P  (spall)  arrival 

a  ^ 

times  and  amplitudes  for  the  explosions. 

Due  to  unknown  complexities  of  both  the  source  environment  and  transmission  path  from 
explosion  to  receiver,  there  has  always  been  a  p''ssibility  that  the  analyses  have  not  really  ex¬ 
tracted  pP  information,  but  rather  some  other  seismic  phenomenon.  In  fact,  many  earthquakes 
show  interference  effects  in  the  first  second  or  two  of  the  P-wavc  train  which  cannot  be  due  to 
a  pP  phase  yet  are  very  similar  to  those  seen  on  explosion  seismograms. 

A  crucial  set  of  data  for  settling  this  questic.t  !.as  been  tabulated  by  D.  L.  Springer.^^  lie 
has  listed  the  pP-P  and  Pg“P  delay  times  for  thirty-seven  U..S.  nuclear  explosions  obtained 
from  accelerograms  recorded  on  the  surface  ove**  the  explosions.  One-way  transit  times  for 
the  explosion  to  the  free  surface  were  essentially  doubled  and  corrected  for  passage  through  the 
fractured  overburden  section  to  obtain  pP-P  times  expected  at  teleseismic  distances.  .Springer 
tabulates  nine  of  the  explosions  which  have  associated  pP-P  delay  times  determined  from  tele¬ 
seismic  seismograms  by  the  methods  described  above.  These  data  are  plotted  in  Fig.  III-ll. 
Also  included  is  an  additional  set  of  seven  pP-P  delay  times  obtained  from  digital  data  at 
NORSAR  by  tliis  author.^ The  ongin  times,  depths,  epicentral  coordinates,  and  available 
yield  information  for  these  explosions  are  given  by  Springer  and  Kinnaman.^^ 

In  Fig.  Ill-il,  the  teleseismic  pP-P  delay  times  for  six  explosions  were  obtained  by  more 
than  one  worker  using  different  data.  Delay  times  for  the  same  explosion  are  connected  by 
vertical  lines  except  for  the  tw'<  points  for  BOXCAR  which  are  bracketed.  With  the  exception 
of  PILE  DRIVER,  most  explosions  in  Fig.  III-ll  show  good  agreement  between  the  surface-zero 
delay  times  and  the  teleseismic  delay  times,  independent  of  the  type  of  analysis  used  on  the 
teleseismic  data.  Generally,  the  teleseismic  delays  are  slightly  larger  than  the  surface-zero 
delay  times.  This  may  be  due  to  the  lack  of  i-esolution  in  time  of  the  teleseismic  data,  or  per¬ 
haps  underestimating  the  decrease  in  P-wave  velocity  of  the  fractured  overburden  section. 

C.  W.  Frasier 
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RELATIVE  AMPLITUDE  (48)  20- SEC- PERIOD  CODA  POWER  LEVEL  (d8) 


Fig.  ni-l(a*b).  Rate  of  decay  of  coda  power  level  and  performance  characteristic 
of  proposed  detection  method  at  ZO-sec  period. 
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Fig.  III-2(a-b).  Rate  of  decay  of  coda  power  level  and  performance  characteristic 
of  proposed  detection  method  at  40*sec  period. 
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Fig.  IIl-5(a-d).  High-resolution  frequency-wavenumber  spectra,  at  40-sec 
period,  for  prefiltered  waveforms  of  27  June  1973  Southern  Sinkiang  and 
Hokkaido,  Japan  events. 
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Fig.  iii”6(a-d).  High-resoiution  frequency-wavenumber  spectra,  at  ZO-sec 
period,  for  prefiltcred  waveforms  of  27  June  1973  Kurile  Islands  event. 
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Fig.  HI-7(a-d).  High-resolution  frequency-wavenumber  spectra,  at  20-sec 
period,  for  prefiltered  waveforms  of  27  June  1973  Southern  Sinkiang  and 
Hokkaido,  Japan  events. 
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Fig.  III-9.  Time  differences  pP-pdP  as  picked  on  VVWSSN  records 
South  American  earthquakes. 
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Fig.  III-IO.  Cross  section  of  seismicity  (mt,  >  5.0)  between  20*S  and  BO’S 
from  Preliminary  Determination  of  Epicenter  lists  for  years  1967-1970. 
Hypocenters  are  projected  onto  a  plane  of  constant  latitude. 


Fig.  111-11.  Teleseismic  pP-P  delay  times 
vs  surface-cero  delay  times  for  nuclear 
explosions  given  by  Springer  (1973). *2 
Each  vertical  line  connects  teleseismic 
delays  obtained  from  different  data  for 
same  explosion. 
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IV.  GENERAL  SEISMOLOGY 


A.  DIGITAL  EVENT  DETECTOR  IN  REAL  TIME 


A  s’ngle-channel  event  detector  \vn'''h  operates  in  real  time  has  been  programmed  for  a 
PDP-7  c(. mputer.  The  purpo.se  of  this  program  is  to  detect  events  that  an  analyst  would  rea¬ 
sonably  pick,  riven  a  single  channel  of  data,  A  direct  application  would  be  the  automatic  detec¬ 
tion  of  events  at  a^.  unmanned  seismological  observar'^ry. 

A  brief  description  of  the  detection  scheme  and  some  examples  with  test  aata  are  presented 
here.  Basically,  the  detection  scheme  compare.*-  .he  spectral  energy  in  e  short  window  with  the 
energy  in  a  long  window  preceding  the  short  window  in  time.  When  the  spectral  ratio  of  short- 
to-long -window  energies  is  greatex'  Lian  an  arbitrary  value,  a  detection  is  declared.  This  is 
similar  to  the  detection  scheme  used  by  SAAC  in  that  short  and  long  data  windows  are  used.^*^ 
Let  ♦he  seismic  tra;.f  Lm  ^  -en  by  the  time  series  x^.  Define  a  moving  periodogram 
at  time  t  given  by 
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where 


At  -  data  sampling  increment  (in  seconds) 
a.’  -  frequency  (in  radian.** /second). 

A  simple  iteration  cap  -e  de.nved  from  Eq.  iIV-1)  to  calculate  the  peric  logranx  at  time  t  +  1, 
namely. 
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From  thesi?  sfjctra,  we  define  short-  and  long-term  average  pc  vers  respectively  by 
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where  I  >  s.  and  L^fu?)  are  averages  of  energies  in  overlapping  windc  /s  of  data.  The 

winoows  used  in  each  a.erage  are  indicated  by  solid  lines  in  Fig.  IV-1.  .-vs  an  event  occurs,  the 
short-term  average  power  incre..se3  before  the  long-term  average  is  affected.  This  suggests 
a  compar'son  between  .sho  t-  and  long-term  aveiages  as  a  detection  criterion.  We  therefore 
calculate  the  average  ratio  of  snurt-  to  long-term  average  for  n  specified  frequencies 


(IV-4) 


A  detection  is  then  declare.^  at  time  t  tf  is  greater  than  a  threshold  pzcameter  f.  If  a 
cietection  occurc,  the  long-term  average  is  fixed  at  its  current  level,  and  the  short-term  average 
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is  continuously  updated.  As  the  event  coda  diminishes,  will  eventually  fall  below  T  and  the 
f-  ector  will  turn  off.  If  no  subsequent  detections  occur  within  i  iterations  past  the  last  detec¬ 
tion,  then  the  long-term  average  will  be  immediately  updated  and  Kq.  {lV-3)  will  again  be  valid. 

We  refer  to  i  as  the  update  lag. 

Experiments  were  run  with  this  scheme  in  order  to  determine  reasonable  values  of  window 
length,  averaging  lengths  and  frequencies  appropriate  for  real  data  from  .NOk.SAH.  In  order  to 
detect  events  with  the  variety  of  periods  observed  at  NOK.'^AH,  four  frequencies  were  chosen: 

1.0,  1.5,  2.0,  and  2.5  Hz.  These  frequenc  are  used  in  the  calculation  of  in  Eq.  (IV-4).  In 
the  examples  presented  here,  the  detection  pai  eters  are  At  -  0,1  sec,  m  =  30,  s  =  30,  and 
t  -  ZOO.  Thus,  the  window  length  for  each  periodogram  is  3  sec,  and  the  short-  and  long-term 
averaging  intervals  are  3  and  30  sec,  respectively.  The  update  lag  i  equals  100,  correspond¬ 
ing  to  a  lag  of  10  sec  after  the  previous  detection  before  the  long-term  average  is  updated. 

To  determine  the  threshold  parameter  T,  the  detection  was  run  on  five  single  channels  of 
noise  at  NORSAR  from  13:02:00  to  13:09:00  fjMT  on  2  3  February  1972.  Values  of  T  =  5.0,  2.5, 
2.0,  and  1.5  were  used,  for  all  five  channels.  The  false-alarni  rate  was  calculated  as  the  ratio 
of  time  the  detector  was  cn  to  the  total  ^ime  processed.  The  combined  results  of  all  five  chan¬ 
nels  are  shown  in  Fig.  IV-2  for  the  fou.-  values  of  T.  The  false-alarm  rate  for  T  =  2.0  is  only 
1.1  percent  but  the  curve  is  rapidly  increasing  that  value.  To  be  more  conservmtive,  the 
value  T  =  2.5  w'as  chosen  for  the  test  cases  desen-ed  below.  The  false-alarm  rate  for  this  T 
is  0.4  percent,  which  is  quite  acceptable. 

In  order  to  test  this  detection  scheme,  beams  of  selected  events  recorded  at  NOR.SAR  were 
added  to  the  five  noise  channels  described  above.  For  each  event,  six  scaled  beams  were  added 
to  the  noise  channels  a  minute  apart.  Each  beam  was  reduced  by  a  factor-of-2  in  amplitude  rela¬ 
tive  to  the  beam  a  minute  earlier,  so  that  a  s^  ries  of  decreasing  beams  was  put  on  each  noisy 
channel. 

Two  examples  are  shown  in  Fig.s.  IV-3  and  IV-4.  In  Fig.  IV-3,  the  event  is  a  shallow  earth¬ 
quake  from  Szechwan  Province  wdiich  occurred  on  16  August  1971.  in  Fig.  lV-4,  the  added  event 
is  a  high-frequency  presumed  explosion  on  6  .Tune  1971  from  Eastei  n  Kazakh.  The  magnitudes 
shown  in  each  figure  are  the  result  of  scaling  the  event  beams  down  by  factors-of-2. 

The  event  detector  was  run  on  the  five  channels  .shown  for  each  event  using  th  ♦hreshold 
parameter  T  =  2,5.  The  results  are  shown  in  Fig.  IV -5,  We  see  that  for  each  evt-.i  the  magni¬ 
tude  4.6  beams  were  detected  on  all  five  channels,  and  two  of  the  magnitude  4.3  beams  of  the 
presumed  explosions  were  detected. 

Figure  IV-6  shows  the  output  of  the  event  detector  for  one  channel  of  earthquake  data. 

Beneath  the  earthquake  channel  is  the  detection  channel  which  Is  either  on  or  off.  The  third 
channel  shows  tne  instantaneous  value  of  the  logarithm  of  R  from  Eq.  (iV-4),  When  R^  rises 
above  the  threshold  valut*  T  =  2.5,  a  detection  occurs  and  stays  on  until  R^  falls  below'  T. 

A  moderate-sized  population  of  events  is  being  studies  to  refine  the  parameters  of  this 

detection  scheme.  r*  c- 

C.  W .  Frasier 

B.  TRAVEL-TIME  ANOMALIES  FOR  THE  GLOBAL  ST.>\TION  NETWORK 

Among  the  products  of  the  aeep-focus  tra\ el-time  .tudy  described  in  .Sec.  ('  below  and  in 
previous  SATS^  are  considerably  improved  estimates  of  the  travel -time  anomalies  caused  by 
variations  ir.  crust  and  upper-mantle  structure  beneath  seismic  stations.  Figures  IV-7  and  IV-8 
show  the  gcog  'aphical  distribution  cf  these  station  anomalies  in  North  America  and  Asii.. 
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One  of  the  main  conclusions  we  have  reached  is  tliat  previous  studies  using  shallow  events 
have  been  contaminated  seriously  by  errors  introduced  by  complex  structures  near  the  sources. 
This  is  show’n,  for  example,  by  a  comparison  of  the  standard  deviations  of  the  data  obtained  in 
the  different  studies.  In  the  largest  previous  study  (Herrin, et  al.'^),  the  standard  deviation  was 
1.00  sec,  whereas  we  obtained  a  value  of  0.f)5  sec.  This  is  in  spite  of  the  fact  that  Herrin,  et  al 
allowed  their  station  corrections  tc.  vary  azimiithally  while  '.ve  did  not.  Indeed,  this  fact  casts 
doubt  on  the  validity  of  the  azimuthal  station  terms  derived  in  previous  studies.  This  doubt  is 
heightened  by  a  comparison  of  such  terms  determined  independently  by  Herrin  and  Taggart,^ 
and  Lilwall  and  Douglas.^  In  both  studies,  the  form  of  the  anomaly  was  taken  to  be 

6t  A  4  H  sin  I?  ‘  01  (lV-5) 


where  f  is  the  azimuth  from  the  station  to  the  event,  and  .4,  H,  and  O  are  constants.  The  dif¬ 
ference  of  two  terms  of  this  form  is  another  function  of  the  same  form,  and  this  fact  ''^n  be  used 
to  examine  the  differences  between  the  two  sets  of  values.  Tor  the  147  stations  com’  s-  to  both 
studies,  the  root -mean-square  magnitudes  (R)  were  0.64  sec  for  Herrin  and  Taggart  and  O.oZ  sec 
for  Lilwall  and  l3ouglas.  The  difference  term,  however,  has  an  rms  value  of  0.68  seel  This 
lack  of  correlation  does  not  mean  that  azimuthal  anomalies  do  not  exist,  but  it  does  in  that 
at  least  most  of  those  previously  determined  are  fjrobably  not  significant  r  the  more 

nearly  error-free  deep-ft>cus  data  can  resolve  azimuthal  station  effects  is  v  being 

studied  „  ... 

R.  K.  .luliai. 

M  Sengupta 


C.  P  AND  S  VELOCITY  STKI  CTI  KE  t)E  THE  LOWER  MANTl.E 

Investigations  of  the  travel  times  of  P-wa.es  from  deef>-focu.s  earthquakes  and  lateral  heter¬ 
ogeneity  of  the  compressional  velocity  distribution  in  the  lower  mantle  have  been  reported  in 
previous  SAT-S.^  These  studies  have  now  been  extended  to  shear  waves  from  twelve  deep-focus 
and  two  intermediate-depth  shocks,  for  which  arrival  time  data  have  been  read  from  70-mm  film 
chips  of  VVVVS.SN  station.^.  The  earthquakes  were  relocated  using  depth  information  from  the 
phases  pP  or  sP  and  the  revised  P  travel-time  curve  and  .station  anomalies.  The  travel-time 
curve  determined  ‘‘or  S-waves  is  shown  i.i  Fig.  lY-^  in  the  form  of  reduced  travel  times  (averaged 
in  Z"  cellst  and  their  standard  errors.  The  solid  line  represents  the  predicted  times  for  a  model 
fitted  to  these  ilata.  Note  that  our  shear-wave  travel-time  data  show  evidence  of  sharp  bends  in 
the  travel -time  curve  at  distances  of  48“,  60 “,  and  71“  which  agree  closely  'with  the  locations  of 
similar  feature.s  found  by  .Johnson^  from  P-wave  dT/dA  data.  However,  since  our  P-wave 
travel-time  data  (which  are  more  accurate)  do  not  show*  these  features,  we  suspect  that  they  are 
fortuitous  results  of  errors.  We  al.so  found  that  station  anomalies  for  .S-waves  are  strongly 
correlated  with  those  for  P- waves,  a  )  are  about  four  times  as  large,  indicating  that  they  are 
due  primarily  ‘o  variations  in  rigiditv  in  the  upper  mantle.  Finally,  in  F’igs.  IV-10  and  IV-11 
we  show  the  regional  anomalies  in  1^  and  .S  travel  times  for  the  distance  range  80”  to  90®. 

Note  the  consistency  between  P  and  S  anomalies  as  well  as  the  fact  that  the  S  anomalies  are 
about  to  6  times  larger  than  the  P's.  Also,  a  striking  feature  in  this  figure  is  the  strong 
gradient  of  the  velocity  anomaly  (especially  for  shear  wave.s)  in  the  neighborhcKid  of  thellaw'aiian 
Islands  for  the  depth  range  of  Z400  to  Z800  km. 

This  work  i.s  presently  being  extended  to  include  amplitudes. 

M.  .Sengupta 
B  ll.  Julian 
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TABLE  IV- 1 

PALEOMAGNETIC  POLES  FOR  RUSSIAN.  SIBERIAN,  AND  EAST  CHINA  PLATESt 

Age 

Russian 

Siberian 

East  China 

Quaternary 

82. 4“  N  146.1“\V 

9  36.2  9.4 

- 

- 

Tertiary 

i 

- 

75.5‘N  38.6*E 

4  630.9  2.8 

Cretaceous 

62.0*N  179. 4*W 

15  27.3  10.3 

7^.9*N  165.3*E 

10  35.7  9.4 

68.0*N  145. 6*E 

*:  35.1  10.9 

Jurassic 

67.0“N  122. 6'E 

4  92.3  7.3 

53.6*N  153.7*E 

4  11.0  22.0 

S3.1*N  151.  3* E 

4  67.9  8.5 

T  riassic 

54.6‘N  153.3“E 

5  100.6  6.4 

41. 5*  N  141.4*  E 

10  17.4  13.6 

- 

Permian 

4b.l*N  167.4*  F 

22  115.7  4.h 

18.2*N  143. 0*E 

9  11.9  16.9 

16.5*N  153.2*\V 

3  19.5  19.6 

Carlwniferous 

42.9“N  170.6*E 

33  16.8  12.8 

31.7*N  145. 4*E 

25  14.0  14.2 

- 

Devonian 

j 

42.4*N  163. 9“E 

6  80.6  6.8 

31.0*N  154. 3’E 

31  33.5  8.9 

- 

Precambrian 

- 

18.0’N  112. 7*E 

7  30.4  10.8 

39.8*N  95.9*E 

4  206.8  4.9 

t  Beneath  each  pole  pos.tion  are  given  the  number  of  data  points  used,  an  estimate  of  the 

Fisher  precision  parameters,  and  the  radius  in  degrees  of  the  95- 
limit  of  the  given  pole,  respectively. 

percent  confidence 
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D.  THK  PAI.KOMACJNKTISM  AND  TKC’ToNlCS  OD  ASIA 


In  an  earlier  SATS  (11  Dec  enil)er  p.  S9),  we  inesenletl  evidenc  e  c)n  the  terlonifs  of 

Asia  from  a  study  of  fault  plane  solutions.  In  that  study,  no  c*leai\  wet  1 -clefined  tectonic  model 
arose  although  seve  1,  non-unicpie  interpretations  in  the  contc^^t  of  plate  tec  tonics  might  he 
suggested  by  the  data  l\ntlier  study  of  fault  |)lane  solutions  and  seisniicitv  patterns  of  centi’at 
and  eastern  Asia  has  temiAed  us  to  suggest  a  tentative  plate  model  for  the  region.  The  b>.  Hil¬ 
aries  of  these  units  the  Kussian.  Siberian,  West  China,  Mast  China,  and  Indian  Plates  -  are 
shown  in  Fig.  tV-12.  Since  paleomugnet ic  data  should  show  evidence  of  these  crustal  units 
moving  independently  in  the  gecdogic  past,  we  have  made  a  study  of  theii-  p'aleornagnetic*  pole 
position  using  published  data. 

The  individual  poles  were  taken  fiom  the  following  authors-  bliramov  and  Sholpo"; 

McElhinny*^;  Irving^®;  Creer^\  VS  ang,  et  tee.  et  al.*  ^  and  1  iu  and  Feng.*"^  f'ctafroin 

rocks  that  were  from  regions  on  or  very  close  to  previously  postulated  plate  boundaries  were 
not  considered  in  the  analysis.  Data  for  which  the  cleaning  procedures  were  not  known  cjr  were 
inadequate  w'ere  also  not  used;  therefore,  rcmsideiable  data  from  the  Kussian  publications  had 
to  be  rejected.  In  general,  only  those  samples  that  were  thermally  cir  magnetically  cleaned 
were  considered.  The  cleaning  methods  used  for  the  Chinese  data  were  unknown,  yet  the  data 
were  used  firstly  because  of  the  paucity  of  that  data  and  seconclly  because  the  ciata  from  widely- 
separated  legions  showed  consistency.  If  there  were  at  least  three  pc^les  for  a  given  age  on  a 
given  plate,  the  mean  pole  and  precision  and  c-ontidenc  e  statistics  were  calculated.  Poles  which 
varied  by  more  than  twice  the  standard  deviation  from  the  mean  pole  were  then  deleted  an»l  the 
mean  pole  was  recalculated  from  the  remaining  data.  Typically,  not  more  than  three  or  four 
for  each  plate  .‘‘or  each  age  had  to  he  rejec  ted  under  this  procedure. 

The  results  of  the  study  are  shown  in  Table  lV-1  where  the  mean  poles  are  tabulated  for  the 
Kussian,  Siberian,  and  East  China  jilates  for  vari»)us  ages.  Paleoniagnetic  data  for  the  Indian 
plate  are  scant,  and  we  found  none  for  the  West  China  plate.  From  Table  IV'-l,  relative  motion 
between  two  plates  would  be  ret  ngnized  by  a  difference  m  mean  pole  position  at  the  same  age. 

It  is  clear  from  Table  IV  -1  that  there  has  been  little  or  no  relative  motion  between  the  Kussian 
and  Siberian  plates  since  Devemian  time  Xlthough  the  I  ral  nu  untains  have  been  suggested  as 
vestiges  of  a  collision  between  the  Kussian  and  Siber  ian  plates,  these  two  units  have  apparently 
acted  as  one  since  Devonian  time  The  Permian  pole  of  the  P^ast  China  plate  differs  from  the 
Siberian  and  Kussian  poles;  however,  the  former  is  based  on  only  three  samples.  Since  Jurassic 
time,  the  East  China  plate  seems  tci  have  expeiienced  little  relative  motion  with  respect  to  both 
the  Siberian  and  Kussian  plates. 

Much  of  the  seismicity  of  central  Asia  e.xists  within  or  near  the  borders  of  the  West  China 
plate.  It  is  under  compression  on  the  southwest  in  the  Himalayan  fold  system  and  under  tension 
in  the  northeast  in  the  Baikal  Kift  Zone  j\dditionally.  from  the  analysis  above,  it  is  bounded 
on  twcj  of  its  three  sides  by  stable  blocks  which  have  not  moved  relative  to  each  other  since  at 
least  the  Jurassic.  If  the  West  (’hina  plate  does  exist  as  a  separate  tectonic  unit,  it  wrould  seem 
to  be,  at  present,  rotating  clccKwise  about  a  pole  somew'here  near  the  southern  tip  of  Lake  Baikal. 

S.  Das 

J  K.  Filsop 
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Fig.  IV-1.  Schematic  of  seismic  data  used  in  detection  process.  Solid  lines 
indicate  overlapping  windows  from  which  each  periodogram  is  computed. 
Windows  indicated  by  S^  and  are  used  in  computing  short-  and  long-term 
spectral  averages  described  in  text. 


Fig.  IV-2.  False-alarm  rate  is  percentage 
of  time  processed  as  a  function  of  detection 
threshold  T.  T  =  2.5  was  chosen  to  process 
test  data. 
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Fig.  IV-3,  Test  array  of  five  noise  channels  at  NORSAR  on  23  February  1971.  starting 
at  *3-02:00  GMT.  Added  to  each  channel  is  beam  of  an  earthquake  from  Szechwan 
Province  1 -minute  intervals,  each  beam  amplitude  scaled  down  by  a  factor-of-2 
from  previous  beam.  Effective  magnitudes  of  scaled  beams  are  indicated. 
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Fig.  rV-4.  Test  array  using  same  noise  channels  as  in  Fig,  Added  to  each  channel 

is  beam  of  presumed  explosion  in  Eastern  Kazakh,  scaled  down  by  a  factor-of-2  every 
minute.  Effective  magnitudes  of  scaled  beams  are  shown. 


16  AUGUST  1971 


6  JUNE  1971 


IsJEChwAN  PRO'GNCE  eastern  KAZAKH 

Fig.  IV -5.  Results  of  procepsing  test  data  of  Figs.  IV -3  and  IV 
detector.  Five  detections  ,.r  a  given  indicate  that  beam 

on  all  channels. 
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P-wave  station  anomalies  in  North  America. 
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Fig.  IV-8.  P-wave  station  anomalies  in  South  Asia. 
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Fig.  IV -9.  Reduced  travel  times  (averaged  in  2'  cells)  of  S-waves 
ind  their  standard  errors.  Solid  line  represents  predicted  t:  -.les 
’or  model  fitted  to  these  data. 
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Fig.lV-10.  Hegional  anomalies  m  P-wavr  travel  times  for  distance  range 
of  80*  to  90'.  .Symbols  are  plotted  at  midpoint  of  great  circle  path.  Crosses 
"epresent  late  arrivals,  and  squares  indicate  crrly  arrivals.  Straight  lines 
are  in  direction  of  source  to  station  azimuth,  their  length  being  proportional 
to  travel-time  anomalies.  Smallest  length  in  this  figure  represents  an  abso¬ 
lute  value  of  0.5  sec,  and  largest  is  an  absolute  value  of  3  sec. 
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Fig.  IV-11.  Regional  anomalies  for  S-wavcs.  Symbols  same  as  in  Fig.  IV-10. 
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V.  COMPUTEH  FACILITIES 


A.  AHPA  NETWORK  COMPUTER  FACMLITV  DEVELOPMENT 

The  ARPA  Nuclear  Monitoring  Research  Office  has  initiated  a  research  and  development 
program  to  make  use  of  satellite  oommuri<*ations,  the  ARPA  Network,  computer  netting,  and 
mass  store  technology  to  achieve  their  seismic  detection  and  discrimination  research  objectives. 
Lincoln  Laboratory  has  undertaken  to  develop  and  implement  special  hardware  and  software 
facilities  which  will  allow  us  to  participate  fully  in  the  program  and  to  exploit  the  full  potential 
of  computer  netting  and  remote  data  bases  for  seismic  research. 

A  survey  of  our  needs  and  technological  alternatives  indicated  that  the  first  step  should 
include: 

(1 )  A  sn;all  host  computer  conne«  ted  to  the  .ARPA  Network  and  having: 

(a)  A  collection  of  basic  peripherals;  card  reader,  printer,  plotter,  mag¬ 
netic  tape  drive,  secondary  storage,  etc. 

(b)  A  number  of  basic  terminals  such  as  teletypes. 

(c)  A  number  of  "intelligent”  terminals,  suitable  for  interactive  graj^ics 
and  typically  containing  a  programmable  minicomputer. 

(d)  A  software  system  prrviding: 

Standard  terminal  access  to  network  resources. 

An  extended  file  transfer  capability  (among  and  between  the 
local  peripherals,  including  the  basic  and  "intelligent”  ter¬ 
minals,  and  remote  machines)  for  high-speed  transfer  of 
data  and  programs. 

Remote  Job  Entry  capability. 

Full  supp<3rt  for  "inJelligent"  terminals. 

(2)  The  re-establi.shment  of  general  Network  capability  on  the  main  Lincoln  Laboratory 
computer  when  the  new  IFaI  370/168  is  installed  in  January  1974. 

(3)  The  creation  of  special  data  access  and  manipulation  software  to  enable  the  seismic 
researcher  and  his  program  to  handle  data  anywhere  in  the  system  including  the 
Datacomputer. 

A  substantial  beginning  has  been  made  toward  these  goals.  The  host  computer,  a  Digital 
Equipment  Corp^iration  PDP- 11/40,  has  been  installed  with  a  card  reader,  an  electrostatic 
printer/plotter,  one  magnetic- tape  drive,  a  1.2- megaword  disk,  and  two  terminals  (an  LA30  and 
a  GT-40).  The  LA 30  is  a  representative  basic  terminal,  capable  of  operating  at  30  characters 
per  second.  The  GT-40  is  a  prototype  "intelligent"  terminal.  It  contains  a  PDP- 11/05  proces¬ 
sor  and,  in  our  configuration,  includes  a  minidisk,  CRT  display  with  light-pen  and  analog-to- 
digita)  conversion  capability. 

A  wide- band  phone  line  has  been  installed  between  the  host  machine  and  the  CC.\  Terminal- 
IMP.  The  Terminal- IMP-to- phone- line  hardware  is  installed;  telephone 'line- to- host  interface 
hardware  is  on  order. 

There  is  a  contract  with  the  University  of  Illinois'  Advanced  Computation  Center  for  a  cus¬ 
tomized  version  of  their  ANTS  Mark  11  system.  This  will  provide  all  the  software  facilities  cited 
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in  (d)  above.  Lincoln  Laboratory  has  been  working  closely  with  the  University  of  Illinois  to 
insure  adequate  support  for  the  "intelligent"  terminals  and  a  sufficiently  powerful  and  flexible 
file  system. 

The  Computer  Systems  Group  at  Lincoln  Laboratory  has  undertaken  ihe  task  of  supplying 
Network  software  for  the  main  Lincoln  machine  (the  370/168).  Initially,  they  will  support  the 
Net  at  the  NCP  and  Telnet  levels.  This  will  give  a  network  user  with  a  basic  terminal  the  same 
power  and  facilities  available  to  the  regular  VM  (time- sharing)  user. 

Our  task  will  comprise  designing  and  implementing  the  special  software  required:  the  Data- 
computer  access  mechanisms;  the  file  transmission  over  the  Net  from  the  370;  etc.  To  make 
this  task  easier,  extensive  explorations  of  SUE  -  a  higher  level,  systems  implementation  lan¬ 
guage  suitable  for  the  370  and  for  the  DEC  machines  (the  11 /40  and  the  GT-40)—  are  under  way. 
In  addition,  bootstrapping  systems  have  been  installed  on  the  two  DEC  machines:  DOS  (Disk 
Operating  System)  on  the  11/40,  and  CAPS  (CAssette  Programming  System)  on  the  GT-40. 

Y.  Peduel 
R.  T.  Lacoss 
J.  R.  Walton 


B.  PDP-7  COMPUTER  SYSTEM 

The  two  PDP-7  systems  continue  to  be  heavily  used  by  Lincoln  Laboratory  staff.  An  appre¬ 
ciable  portion  of  the  remaining  time  is  used  by  visiting  seismologists,  by  graduate  students,  and 
for  ISM  production  runs. 

The  system  has  become  fairly  stable,  and  it  was  decided  ihat  a  complete  revision  of  system 
documentation  was  in  order.  The  new  Software  Handbook  will  consist  of  twelve  chapters  in  two 
volumes.  Of  these,  nine  chapters  have  been  written  and  published: 

I.  System  Introduction 

II.  Basic  System  Utilities 

III.  Editor 

V.  Assembler 

VIII.  Data  Analysis  Console 

IX.  Data  Analysis  and  Display  System 

X.  Program  Library 

XI.  Subroutine  Library 

XII.  Preliminary  Index 

Yet  to  be  published  are: 

IV.  Fortran 

VI.  I/O  Devices 

VII.  Data  Formats 

L.  J.  Turek 

M.  F.  O'Brien 
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